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(57) Abr6ge/Abstract: 

The operation of a fossil-fueled thermal system (20) is quantified by obtaining a reference fuel chemistry before on-line operation 
(120). and thereafter operating on-line. In on-line operation (250). a set of measurable operating parameters is measured, 
including at least effluent concentrations of oxygen and carbon dioxide, and optionally the concentration of effluent water and 
the concentration of effluent sulfur dioxide. An indicated /Vir/Fuel ratio is obtained, as are the ambient concentration of oxygen, 
and air pre-heater (36) leakage (40), and dilution factors. The fuel ash (Eq.(54)) and fuel water (Eq.(47)) are calculated, and the 
complete As-Fired fuel chemistry is calculated. From the complete As-Fired fuel chemistry (Eq.(29)), the pertinent systems 
parameters such as reference fuel heating value (Eq.(65)). boiler efficiency (Eq.(18)), system efficiency (Eq-(3)), fuel flow rate 
(Eq.(21)), total effluent flow rate (Eq.(71)), individual effluent flow rates (Eq.(69)), and individual emission rates (Eq.(70B)) are 
determined in a fully consistent manner. 
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INPUT/LOSS METHOD FOR DETERMINING 
FUEL FLOW. CHEMISTRY. HEATING VALUE AND PERFORMANCE 
OF A FOSSn>-FIRED SYSTEM 

ABSTRACT OF THE DISCLOSURE 

The operation of a fossil-fiieled thermal system (20) is quantified by 
obtaining a reference fuel chemistry before on-line operation ( 1 20), and thereafter 
operating on-line. In on-line operation (250), a set of measurable operating 
parameters is measured, including at least effluent concentrations of oxygen and 
carbon dioxide, and optionally the concentration of effluent water and the 
concentration of effluent sulfur dioxide. An indicated Air/Fuel ratio is obtained, 
as are the ambient concentration of oxygen, and air pre-heater (36) leakage (40), 
and dilution factors. The fuel ash (Eq.(54)) and fuel water (Eq.(47)) are 
calculated, and the complete As-Fired fuel chemistry is calculated. From the 
complete As-Fired fuel chonaistiy (Eq.(29)), the pertinent systems parameters 
such as reference fiiel heating value (Eq.(65)), boiler efficiency (Eq.( 1 8)), system 
efficiency (Eq.(3)), fiiel flow rate (Eq.(2l)), total effluent' flow rate (Eq.(71)), 
individual effluent flow rates (Eq.(69)), and individual emission rates (Eq.(70B)) 
are d^rmined m a fiilly omsistent manner. 
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OPERATION OF A FOSSE^FUELED THERMAL SYSTEM 



BACKGROUND OF THE INVENTION 

5 This invmtion relates to a fossil-fired power or steam generation thermal 

system, and, more particularly, to a method for determining its fiiel chemistry, fiiel 
heating value, fuel flow, and thermal performance torn its basic operating 
parameters. 

The importance of accurately determining thennal efficiency is critical to 

1 0 any thermal system. If practical day-to-day improvements in efficiency are to be 
made, and/or problems in thamally degraded equipment are to be found and 
corrected, then accuracy in determining thermal efficiency is a necessity. The 
tracking of the efficioicy of any thermal system lies fundamentally in measuring 
the useful output, and the total in-flow of fiiel. The useful output from a fossil 

15 fueled system includes the generation of electrical and/or mechanical power 
and/or the production of a heated working fluid such as steam. 

The measuring of the usefiil output of thmnal systems is hig^y developed 
and involves the dn^ measuremmt of electrical output and/or mechanical drivels 
and/or thermal energy flow. Measuring electrical and mechanical power is well 

20 established. Measuring th^mal energy flow involves direct measurement of the 
inlet and outlet pressure, temperature and/or cpiality of a fluid being heated by tfie 
combustion gases, as well as measurement of its mass flow rate (m). From this 
information enthalpies (h) may be determined, and thus the total energy flow, 
ni(h^ci • delivered from the combustion gases, wh«i also considering 

25 incidental losses, may be determined. 
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The measurement of the energy flow of the input fiiel requires knowledge 
of the heating value of the fuel and its mass flow rate. For thermal systrais using 
gaseous fuels, tiie fuel's composition may be weU characterized, thus its heating 
value may be determined based on known heats of combustion associated with 
S individual components. 

However, there are numoous situations v^ere a fossil-fiieled system's fuel 
eneigy flow is not well characterized. For example, even a gas-fired system, 
having no on-site fuel gas analyzer, may receive fuel from multiple sources whose 
composite heating value variation is ±5 percent or greater. The measurement of 

10 fiiel flow may often present a problem of measurement accuracy, especially at 
lower loads v/here flow measurement devices are not calibrated. In either case the 
determination of fuel energy flow is critical for proper thermal understanding of 
any fossil-fired syston, either for direct confumation of computed results and/or 
for improving system thermal efficiency. 

IS The importance of accurately determining pollutant concentrations and 

their flow rates is also oitical to the practical operation of any fossil-fired system 
due to en^dronmental constraints imposed through regulation, the potential of 
regulatory induced fines, and concern by the owner of the facility for 
en^dronmental protectioiL 

20 Given these considerations, it is equally important to have analytical 

evidence of the enxirs made in the detennination of fiiel flow and the errors in fiiel 
heating value, and thus the resultant errors in the determined thermal efficient 
of the system. Further, any error in effluent flow, which is proportional to any 
error made in fuel flow, is significant when monitoring and reporting to regulatoiy 

25 agencies the effluents firom any source of fossil combustion. The method of this 
invention provides a tedinique for specifying and correcting possible errors based 
on the consistency of the primary measurements of effluent O^, COj, indicated 
Air/Fuel ratio, assumed or measured effluent HjO, possible air pre-heater leakage 
and the concentration of O2 in combustion air 

30 The measurement of fuel flow has traditionally been accomplished via 

measurement of its mechanical effects on a device. Such effects include the 
pressure drop across nozzles or orifice plates, unique fluid densities, integrated 
weighing of a fiiel handling conveyor belt (commonly used for bulk fiiels such as 
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coal), speed of sound, nuclear resonance, chaise in bulk storage levels, etc. 
Present industrial tediniques for measuring ps or oil fuel flow result in typical 
errors of 1 percent to 10 percent relative to true values, depending on the care 
taken in designing, manufacturing, installing and calibrating the flow metering 
5 equipment and in its data reductioa Under ideal circumstances, titter accuracies 
(i.e., smaller errors) are possible for gas and oil fuels, reaching at best 0.2S 
percent, but diis is consid^ed very unusual, always requiring ^ctraordinaiy 
e^qsense. 

For bulk fuel such as coal, bio-mass, slurry fuels, wood, agricultural 

10 byproducts such as shells 6rom nuts, trash and refuse, the typical accuracies of 
flow metering range upwardly fiom 5 percent and higher. Historically, bulk flow 
measurements have such poor accuracy that they are used only as a relative 
indicator of fuel flow. For fossil-fired systems any fiiel flow error greater than 
appro?umately 1 percent, and certainly greater than 2 percent, is sufficiently high 

IS to preclude trending of the monitored fuel flow rate for reasons of thermal 
efificiency or for detecting degraded equipment Improvement of efficiencies in 
a thermal system is classically concemed with a number of small incremental 
improvements, typically each in the range of 02 percent to 0.6 percent A dozen 
or more of these, taken together, may result in 3 percent to 6 percent 

20 improvement For example, an average 4.5 percent improvement has been 
physically demonstrated at over two dozen conventional power plants, see 
Lang, "Methodology for Testing and Evaluating Poww Plants Usmg Computer 
Simulators", 1990 ?edomaTice Sofbvare User's Groiq> Meeting, May 1-4, 1990, 
St Louis, sponsored by EI International, now Scientech Inc. of Idaho Falls, ID. 

25 Prior ifiproaches ^ch attempt to address the accurate determination of fuel flow 
are discussed below. 

Another critical consideration in determining thermal efQciency is the 
variation in the fuel's heating value due to variations in fuel chemistry. Chemical 
variations appear through the mix of fuel water, fuel mineral matter (called fuel 

30 ash), and the relationships of the elements comprising the basic hydrocarbon 
molecular chain and any free inorganic elemmts: nitrogen, oxygen, carbon, 
hydrogen and sulfur, but principally carbon, hydrogen and oxygen. If an accurate 
and direct flow measurement of bulk fuels is not practical, the only alternative is 
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the detennination of fuel energy flow, which is the product of flow rate and 
heating value of the fuel, based in part on the measured energy flow to the 
working fluid. If errors exist in the heating value, either an assumed, measured, 
or calculated value, errors will then result in the fuel flow. Prior approaches 
5 which attempt to address the determination of fiiel energy flow are discussed 
below. Further, over the past 140 years of producing safe high pressure steam for 
society, mostly fix>m coal, there has been no invration or process even suggesting 
an ability to determine a coal's fiiel ash content in real time based on 
thermodynamics. 

1 0 The ^proach of this invention is a much imiMOved "Input/Loss Method'*. 

Prior input/loss methods have been known to the inventor, and to T. Buna as early 
as 1955. 

One prior approach related to the present invention was developed by T. 
Buna m 1 955 for the analysis of multiple fuels fed to a power plant His approach 

15 was to characteri2se a fuel's effluent CO2, given differing effluent O2 values, by 
assuming fuel diemistiy of the individual fiiels. With tfiis data for multiple fuels 
and knowledge of ttie Useful Energy Flow Delivered, he advocated determining 
eadi fuePs flow rate. He presents an ''output-loss"* and an "input-loss'' approaches 
to determining boiler efficiency. The present invoition is related only in that a 

20 course reversed from Buna's method is accomplished. This invention computes 
fiiel diemistiy based on effluent measurements, it assumes that all secondary 
fuels, unlike Buna, are known having defined dittiistries, heating vahies and As- 
Fired flows. The reference is: T. Buna, "Combustion Calculations for Multiple 
Fuels*', ASME Diamond Jubilee Aimual Meedng, Chicago, IL, November 13-18, 

25 1955, Paper 55-A-185. 

Another related art to the present invention was developed by the Electric 
Power Research Institute (EPRI) at the Morgantown power plant, a coal-fired 
conventional system. This technique, termed the "Output/Loss" Method, is 
described by E. Levy, N. Sanmac, H.G. Grim, R. Leyse and J. Lament, 

30 "Output/Loss: A New Method for Measuring Unit Heat Rate", Am. Society of 
Mech. Engrs., 87-JPGC-Pwr-39. The Output/Loss Method produces boiler 
efficiency {r\^ indepaidoit of fuel flow. Assuming a conventional power plant, 
by determining the mcrgy flow to the working fluid (Em Ah) and gross electrical 
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power production (P)» systen thermal efficiency may be determined, i.e., 
T)teiie/'/(^^^h)' addition, atttiough it is not the objective of the Output/Loss 
Method, if heating value and boiler ena:gy credits (HHVP + HBC) are known, 
then the As-Fired fiiel flow (m^) may be determined as m^ = 

5 {2jnAh)/[iiboaa(HHVP + HBC)]. Use of boiler energy credits (HBC), and gross 
power (P) versus net power, are discussed in the Preferred Embodiment. The 
technique relies on measuring emission gas flow directly, and thus 
Knowing emission gas flow allows the determination of the majority of the 
thermal losses associated with combustion, called ''Stack Losses". However, this 

0 approach has drawbacks when it is applied for practical applications in power 
plants for the following reasons: 1) it does not rely on measurement of flue gas 
concentrations (thus changes in the in-flows of water/steam to the combustion 
process, or changes in effluent CO2, as might reflect changes in fuel chemistry or 
air pre-heater leakages); 2) the errors in effluent gas flow measurements in 

S irregular ducts not designed for accurate flow measurements, vdiich is the case at 
most power plant fiacilities, may easily exceed i20 pox^ent, resulting m over ±4 
percent error in systm efficiracy since typically stack losses are «20 percoit (i.e., 
±4 percent error in fuel flow); 3) the technique of direct flue gas flow 
measurements does not consistently meet current U.S. &ivtronmental Protection 

0 Agency's accuracy requirements of ±15 percent; 4) it is obvious from the above 
discussion that if the fuel's heating value (HHV) is variable, as is common with 
most coal-fired systems, and is not properly monitwed in a continuous manner, 
then calculated flxel flow will also be in mor, due to mors made in the assumed 
heating value; S) direct measurement of effluent flow commonly involves 

S ultrasonic, sonic or direct gas velocity measurements, requiring assumpti<ms as to 
effluent compositions, i.e., fuel constituents and air pre-heater leakage, errors in 
these will force mors in the deduced effluent flow; and 6) if the fuel bears highly 
variable quantities of ash (a pure dilutive or concentrative influence on fuel 
heating value, but affects effluoit flows through Air/Fuel relationships), then the 

0 computed fuel flow, since it is dependent on an assumed heating value with an 
assumed ash content, will likely be in error. 

Another related art to the preset invention was presented at a technical 
conference in 1988 by S.S. Munukutla, et al.. In this work the authors develop a 
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process which used effluent measurements to determine tfie Moisture-As-Free 
(MAF) composition of coal, their published work teadies little relative to the art 
of monitoring theimal p^onnance and the determination of coal chemistiy of 
practical operating powa- plants. Munukutla, et al, do not consider air pre-heater 

S leakage, but measure combustion efQuents before the air pre-heater (at the 
economizer outlet). They invoke high accuracy efiQuent measurements afforded 
with gas diromatography, not common industrial instruments. This basic 
technique, commonly referred to as Thermal Analysis using in part gas 
chromatography, is used in laboratories to det^mine heating values under highly 

0 controlled conditions, refo: to Chapter 9, **Measurement of Heat of Combustion" 
contained in Steam. Its Generation and Use. 40th Edition, edited by S.C. Stultz 
and J.B. Kitto, published 1992 by the Babcock & Wilcox Company, Barberton, 
OH. It was also the subject of inventor's earlier US Patent S3273S6. They do not 
consider the injecti<m of limestone as is common with fluidized bed combustors, 

5 representing the largest single type of steam generator sold out^de the U.S. They 
make no claims for the on-line detenninaticm of fuel ash, but assume it is constant 
based on laboratory analj^is. Eirors are made in their Equation 3, which if 
applied to their process would yield ridiculous results. This error is corrected in 
subsequent works by Munukutla as presented in 1989 and later. The ratio of 

0 atmospheric non-Oj (principally N2andAr)to02issetasa constant in their woric 
at 3.76 (vs4iich assumes O2 at 21.0 percoit). First, this ratio, commonly found in 
fundamental thermodynamic text books, is in error, the coirect ratio is 3.7737, 
herein termed (j)^^^ derived from a value of 20.948 percent atmospheric Oj per 
NASA (reference U.S, Standard Atmosphere 1976,NOAA-S/T-76-1562-NASA), 

5 Not using a correct value results in an error in determining excess air and molar 
oxygen balances. Second, in iq}proximately 1 0% of the power plants tested by the 
inventor, atmospheric O2 may be degraded, resulting in a higher 4>Act value, due 
to in-leakage of flue gases into the combustion air stream and/or local atmospheric 
inversions which may result in starving the local oivironment for oxygen. In the 

} present invention this ratio is expressed as a variable, (}>^ to be set by the user 
based on circumstances local to the thermal system. 

At a systCTis un<testanding level, the 1988 approach by Munukutla, et al., 
is lacking for they compute boilra* efiBciency only as a final result of the technique 
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using the Input-Output Method (discussed below); see their Equation 10. A 
computed boiler efficienQr incorporated within the technique proper is lacking, 
thus integral consistency reqmred f<^, and a feature of, the present invention is not 
assured using their methods. A most important attribute of the present invention 
S demonstrates an integral relationship between coal constituents and the 
determination of boiler efficiency, leading to system efficiency; indeed without 
coal constituents boiler efHciency simply cannot be computed by the present 
invention. Munukutla, et al.^ technique determines tfie total fuel enagy flow from 
system energy flow balances (their Fig, 2 and Equation 4), developing a "fuel 

10 energy" term, Am^^p = + Qst«mi; offering no explanation as to how the term 
"A" in this e?q>ression is related to coal heating value, as one might assume. 
However, this expression is then used to develop mass flow rates of the coal*s 
constituents, termed fij (in pounds/second). These terms are tiien used in a 
correlation which develops total fuel meigy flow (Btu/hr), thus bypassing a boiler 

IS efficiency calculation. Air pre-heater leakage is not considmd in their system 
energy balances. It is noteworthy that this energy balance technique, leading to 
their Equation 4, is not referenced in later publications on the same subject by tfie 
same principal author (in 1989, 1991 and 1995), however the basic 
stoichiometrics, corrected, are referenced. The reference for this work is S.S. 

20 Munukufla, G. Tsatsaronis, Y. Shih, D.E. Anderson and S.M. Wilson, "A 
Microcomputer Software for On-Line Evaluation of Heat Rate", ASME Power 
Generation Conference, Philadelphia, PA, September 25-29, 1988, Paper 88- 
JPGC/PTC-L 

In three later works, Munukutla and his colleagues attempt to refine the 
25 technique. At a 1989 conference, published in 1990, Munukutla and A. Bose 
again propose using gas chromatography to measure economizer outlet gaseous 
effluents in determining MAP heating values. The heating value correlation 
employed, using weight fractions of coal's constituents, was the Mott-Spooner 
correlation based on Dulong's formula. These are both well known correlations 
30 in the industry, however they are not based on chemical binding enopes as this 
invention employs. However in this 1989 wock« the correlation was used 
incorrectiy; the correlation was intended for MAP fuel constituents, Munukufla 
and Bose inconsistendy used the wet As-Fired. As witness to the inaccuracy of 
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their work, their calculated MAF heating value reasonably agreed with the As- 
Fired (wet). In 1991, Munukutia, P. Chodavarapu and D.C. O'Connor published 
essentially the same woik as in 1989, again based on gas chromatography effluent 
measurements. Intheir 1991 work fuel wateriseithermeasuredordeterminedby 
S differmce, assuming the remaiiung effluents are measured on a wet base. It is 
obvious fincHOi their work that reasonably accurate molecular weights are not used; 
given the sensitivity of any such approach, such enws amount to 1% error 
(using their numerical examples, the fuel carbon fraction is computed as 0.6616, 
versus the conect 0,6570 as determined by simply using correct molecular 

10 weights. Further, as with the 1988 and 1989 work, no provisions are made for 
limestone injection. Again the 1991 work used the Mott-Spooner correlation 
based on Dulong's formula to determine the heating value, HHVmaf> but this time 
the authors used the consistent MAF fuel constituents, later correcting HHV^^ 
for fuel water and ash. The result was in error with the reported value by 387 

1 5 ABtu/lbm, or 3. 1%. Using only their 1 99 1 reported fuel analysis, this invention's 
methods produced an error of 162 ABtu/lbm. They also reported results of 14 
tests preformed, resulting in a standard deviation of ±785 ABtu/lbm, or ±5.6%. 
In 1995,MunukutlaandF. Khodabakhsh published a shnilar work as in 1 989 and 
1991, but this time based effluent measurements on continuous emission 

20 monitoring system (CEMS) instnmientation. CEMS instrumentation is required 
by the US EPA on atl stationary sources of pollutants. Effluent measurements 
include CO2, CO, Oj, SO, and volumetric flow rate. Undsr US and European 
regulations these measurements must be made at the boundaiy of the system, at 
the air pre-heater's outiet, i.e., the smiE^e stack, not at the boiler's outiet (i.e., the 

25 air pre-heater inlet); refer to FIG. 6. Munukutla and Khodabakhsh s^plied the 
same stoichiometrics for this work as was used since 1 988, applicable for a system 
without air pre-heater leakage. Given the measurements are made with air leakage 
present, die method is fundamentally flawed; all effluents emanating from the 
combustion process are of course diluted with air leakage. However, the system 

30 was defined in such a manner that air leakage was factored out when considering 
a system energy balances. This was done by using the air pre-heatefs hot outlet 
flow to the boiler, per their Figure 2, not the combustion air inlet to the air pre- 
heater. The cwisequence of this is tiiat although fuel flow could be computed 
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correctly, the integral det^ination of boiler efiBcienqr, as with the 1988 work, 
is not possible. Further^ since the proposed effluent measuranents in this 1995 
work are inconsistmt with its basic analytics, the computed fuel constituents are 
then seriously flawed Limestone injection is not considCTed. Fuel ash is assumed 
5 constant In all of these works, from 1988 through 1995, no mention is made of 
converting the resultant heating value from a constant volume base to a constant 
pressure base as requued given that laboratory determined coal heating values 
derive from "bomb calorimeter** devices, a constant volume process. However, 
a refinwnent described in their 1995 work suggests an iterative correction of the 

10 measured CEMS volumetric flow with that computed through stoichiometrics. 
The present invention does not employ effluent flow measurements and thus this 
reflnement is not mat^al. 

RefM-cnces for the 1989, 1991 and 1995 woiks by Munukutla and his 
colleagues include: S. Munukuda and A. Bose, '*Qn-Line Elemental Analysis of 

15 Coal Using Gas Chromatogrq)hy", 1989 EPRI Heat Rate Improvement 
Conference, May 3-5, Knoxville, TN, published 1990); S. Munukutla, P. 
Chodavarapu and D.C O'Connor, **Qn-Line Coal Analysis from Measurement of 
Flue Gas Components", ASME International Power Generation Conference, San 
Diego, CA, October 6-10, 1991, Paper 91-JPGC-Pwr-17; and S. Munukutla and 

20 R Khodabakhsh, ^'^ihancraient of Boiler Poformance Evaluation Methods Using 
CEMS Data**, ASME International Joint Power Generation Conference, 
Minneapolis, MN, Oct 8-12, 1995, PWR- Volume 29, 1995. 

Another approach was developed by the present inventor and was 
described in United States Patent 5,367,470. The invention of the *470 patent is 

25 noteworthy because it addresses the detennination of boiler efficiency without 
knowledge of fiiel flow and without knowledge of effluent flow, but knowing 
principally effluent CQj, effluent HjO and routine system data, and it is usefid in 
many situations- However, the Air/Fuel ratio is not reqiured in the metiiod of the 
*470 patent. A key to the invention of the *470 patent is its requirement to 

30 repetitiously adjust, or iterate, on "an assumed water concentration in the fiiel until 
consistracy is obtained between the measured CO2 and H2O effluents and those 
determined by stoichiometrics based on the chemical concentration of the fuel". 
Some aspects of the invention are dependent upon high accuracy effluent water 
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and carbon dioxide concentration measurements, or "predeteimined accuracies" 
of these effluents. The present invration has no sudi limitations on accuracy. The 
difficulty in the method of die '470 patent lies with the fact that adjusting fuel 
water, whidi will of course alter the computed effluent water, has no prima facie 
5 effect on a dry-base effluent CO^. It is true, for exanq^le, that if fuel water is 
increased, the relative fraction of the other fuel's constitu^ts, per unity mole of 
total As-Fired fuel, will decrease assuming that the fuePs other constituents, 
nitrogen, oxygen, carbon, hydrogen, sulfur and ash, remain proponionally 
constant to each other. However, it would be unusual diat any given fuel water 

10 adjustmOTt would produce an exactly consistent effluent COj, with the exception 
where the dry chemistry is constant or highly predictable. Further, if the fuel has 
a variable ash content, ash having a pure dilutive or concentrative influence on 
fuel chemistry and fiiel heating value, then such variable effect coidd not possibly 
be determined by merely iterating on fiiel WBto-. A higher assumed fuel water 

IS may decrease a wet-base effluent COj, but the actual fuel could contain much 
lower ash, thus actually increasing die amount of fuel carbon relative to the v^ole. 
The approach of simple water iterations of die '470 patent is useful in many 
situations, such as ^ere the coal fuel bears littie and constant ash, and, further, 
where v^ high accuracy and consistmt effluent CO2 and H2O measurements are 

20 made, but has iinutations in otiier applications. 

Yet another approach was developed by the present inventor and was 
desCTibed in United States Patent 5,790,420. The invention of the '420 patent is 
noteworthy as it extends the sq^proach of the '470 patent to include combustion 
turbine systems. The '420 patent is concerned with mediods for improving 

25 themial efficioicy, det^mining effluent flows and determining fuel flow of fossil- 
fired systems through an understanding of the total energy flow (fitel flow rate x 
heating value), the As-Fired input to the system. The '420 patent states that 
method errors will offset "the sign of the error introduced by the heating value 
will always have an opposite change in the calculated fuel flow". Errors may be 

30 introduced in the '420 patent by the use of its Equations 3 1 or 32 to compute the 
dry-base heating value, dependent on knowledge of the dry molar coii^>osition of 
the fuel. These dry compositions niay be determined throiigh "use of a correlation 
relating carbon, hydrogen, oxygra and sulfur contoits to a diy-base heating value 
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film conrectijig for water". Fuel ash is suggested by the *470 patent and the '420 
patent as being treated as a constant value. This is the case since the fuel 
constituents are solely defined per unity moles of diy-base fuel or as an As-Fired 
(for wet-base) fuel. A diy-base fuel contains ash; a wet-base fuel contains ash and 
5 water. Ifthe effects ofvariable ash were to be addressed, fuel constituents would 
by necessity initiate fix>m a base free of both ash and water, i.e., so-called 
"Moisture-Ash-Free"; such a base is not mailioned in eitho- Ae *470 patait or the 
*420 patent The *420 patent explains that the molar quantity of fuel water "is 
iterated until convergence is achieved", resulting in an As-Fired heating value. 

1 0 Again, as water is altered, the aggregate of all other fuel constituents are altered 
in opposite fashion to maintain a normalized unity moles of fuel. As with the 
approach of the *470 patent, the *420 patent requires high accuracy 
instrumentation, stating "the apparatus necessary for practicing the present inven- 
tion includes utilization of any measurement device which may determine the 

IS efiQuent concentrations ofHp and CO2 to high accuracr/**. The i^roaches of the 
'420 patent and the '470 patent, which are rudimentary Irq^ul/Loss methods, are 
dependent on thermodynamically understanding a fossil-fired system without 
direct measurement of fuel or effluent flows, but there is room for improvement 
Another related art to the present invention involves using fuel water and 

20 fiiel ash instruments intended for on-line operation. By fuel ash is meant the fuel's 
non-combustible mineral content, before firing. These instruments employ a 
variety of techniques- Fuel water instruments include: capacitance techniques, 
microwave techniques, ultrasonics techniques and IR spectroscopy. Fuel ash 
instruments include: X-ray backscatter. X-ray fluorescence, gamma-ray 

25 backscatter, dual energy gamma-ray transmission, garruna-ray pair production, 
natural gamma radiation, prompt gamma neutron activation analysis, laser 
spectroscope, electron spin resonance and nuclear magnetic resonance. Knowing 
the fuel's water and ash contents can lead to adjustment of the assumed heating 
value. This provided the fuel's MAF chemistry remains constant Further, these 

30 methods do not involve techniques in which errors made in the fuel ash 
measurement are ofif-set by fuel water concentrations. Further, common industrial 
accurades are no better than ±5%, thus a ±5% enor in fuel heating value* Further, 
whereas the direct determination of fuel water and ash would aid the present 
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invention as overchecks of its calculated values, such instruments by themselves 
do not provide an integrated ^proach to the understanding of thermal systems. 
Clearly the computation of boiler efficiency, as integrally related to system 
parameters, is not made by simple use of these instruments. 
5 Another related art to the preset invention is Thmnal Analysis whidi 

employs laboratory techniques of dififermtial thermal analysis (DTA) and 
di£ferential thermal gravimetrics (DTG), also termed thermogravimetry, combined 
with gas diromatQgraphy or other high accuracy gas analyzer, to determine the 
elementary analysis and proximate analysis of coal. This technique is intended for 

1 0 a laboratory envirorunent, since it analyzes only gram amounts of coal . Obviously 
using tons of coal/hour by a thermal system such as a 100 Mwe power plant, 
developing a representative gram-size sample would present an insurmountable 
problem. One of the largest companies which supplies such equipment is Mettler 
Toledo, Hightstowns, NJ which has been selling such equipment since at least 

1 5 1 980. TTiis basic technique is mentioned in Chapter 9, ''Measurement of Heat of 
Combustion" contained in Steam. Its Generation and Use, cited above. 

Another approach was developed by the American Society of Mechanical 
Engineers (ASME) and published as its Power or Performance Test Codes (PTC). 
Several of these codes discuss two methods which are relevant: the Ii4)ut-0utput 

20 Method invoked in PTC 4. 1 (Steam Generators), 4.4 (Gas Turbine Heat Recovery 
Steam Generators), and PTC 22 (Gas Turbines); and the Heat-Loss Method 
invoked in PTC 4. 1 and PTC 4.4. Tlie Input-Output Method relies on the direct 
measuiement of fuel flow. For coal-fired plants or any bulk-fiielqrstems, this has 
no applicability for the improvement of thermal efficiency for tfie above-discussed 

25 reasons of inaccuracy in fuel flow measurements. The Heat-Loss Method was 
intended to address the issue of inaccurate coal flow measurements by 
determination of unique stack and non-stack losses, thus: ^symm = LO - 
E(System Losses). The difficulty with this method, as in the case of the 
Output/Lossmetlu)d, lies in the need tomaketwocriticald^enni^ l)eith^ 

30 measurement of gffseous effluent flow directly or accurately knowing fuel 
chemistry leading to an effluent-to-fuel flow determination; and 2) v/hcn applied 
to coal-fired systems, measuring effluent ash flow and its associated unbumed 
carbon in the ash. As discussed, the required accuracy is not possible on a 
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consistent basis vAisn measuring efiOuent flows. Further, the traditional use of tfie 
Heat-Loss Method requires an "accurate, simple, and ultimate analysis of the fuel 
being fired,** The difficulty in determining both "bottom" and "fly" ash flows, 
defined in PTC 4.1 , is evidoiced by the fact the procedure is not prefomied on any 
5 routine basis by any known coai-fired power plant 

In summary, the ^roaches of the Output/Loss Method, Munukutla and 
his colleagues, the *470 patent, the '420 patent, fiiel water and fiiel ash 
instrumentation, laboratory Themial Analysis, the Input-Ou^ut Method, and the 
Heat-Loss Mediod all have significant limitations. The Output/Loss Method and 

1 0 both ASME Methods are flawed conceptually for at least typical large coal-fired 
systems. None of these methods are ^plicable for on-line monitoring. The 
methods of the *470 patent and the '420 patent do not consider: 1) the 
determination of fuel ash (intrinsically assuming a constant or known relationship 
between fuel carbon and fiiel ash); and 2) the complexities of the non-wat^ 

IS constituents, iterating sunply on fuel water. That is, the methods of the '470 
patent and the '420 patent sinq>ly alter fuel water under the assumption that the 
relationships betweoi the fuePs non-water constituents remain as constants or are 
predictable throu^ correlations. If die determination of the fuel's constituents is 
flawed, then the determination of the fuel's heating value is flawed, and thus the 

20 fiiel flow vdll be in error. Although as the '470 patent and the '420 patent state, 
such errors in heating value and fuel flow tend to ofif-set one another, this is 
typically true only if the ash is both a relatively small fisction of the total fuel and 
of essentially constant concentration. As will be seen in relation to the present 
invention, tolerating such errors voids useful information associated with absolute 

25 knowledge of heating value, voids the accurate determination of effluent flow, and 
voids any computational overcheck of the accuracies of the effluent 
measurements. None of these methods, except the '470 Patent and the '420 
patent, consider multiple fuels which are commonly used in commercial power 
plants. None of these methods considers the thermodynamic determination of 

30 fuel ash. None of these methods considers the use of limestone injected into the 
combustion process. None of th^ methods, except the '470 Patent and the '420 
patent, consido* air pre-heater leakage such that gas concentrations on either side 
of the air pre*heater are useable in stoichiometric relationships. None of these 
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methods, except the '470 Patoit and the '420 patent, consider variable Q2 in the 
combustion air local to the system. 

Another approach was developed by the United States Environmental 
Protection Agency (EPA) as related to the determinations of efHuent flow and 
5 individual emission rates ObpoHutan/mllion-BtUf^. EPA's approach is described 
in its regulations, Chapter 40 of the Code of Federal Regulations (40 CFR). 
Specifics are described in 40 CFR Part 60, Appendix A, Methods 1, 2 and 17 
defining various teduiiques for measuring effluent flows, and in 40 CFR Part 60, 
Appendix A, Method 19 defining "F Factors" used to determine emission rates. 

10 The EPA requires the direct measurement of effluent flow &om stationary sources 
of fossil combustion. The EPA also requires the reporting of emission rates for 
the major pollutants, determined using the F Factor technique. The EPA's 
approach has the same shortcomings as discussed for the Output/Loss Method. 
The EPA q)proach does not require any inter-relationship between a computed 

1 5 fuel flow, which through stoichiometrics must be consistent with effluent flow. 
In a work by Lang, et al.^ reporting actual test results on a large power plant, 
EPA's Methods 1, 2 and 17 resulted in effluent flows, on avmge, 12 percent 
higher than those consistmt with systrai efficiency; this implies a 12% high^ fuel 
flow, a large error in system understanding. Further, all EPA methods produced 

20 higher flows relative to those consistent with system efficiency^ and none were 
consistent in themselves. These results are typical of the 12 to 15% bias reported 
throughout the industry. See F.D. Lang, et al, "Confirmatory Testing of the 
Emissions Spectral Radiomrter/Fuel Flow (ESR/FF) Instrument", Electric Power 
Research Institute (EPRI) 1994 Heat Rate Improvement Conference, May 3-5, 

25 Baltimore, MD. In another study at two large power plants, use of EPA's Method 
2 produced 9.8 percent and 18.6 percent higher system heat rates (i.e., system 
efficiencies) based on measured effluent flows, see R.D. McRanie, et al, "The 
Electric Power Research Institute Continuous Emissions Monitoring Heat Rate 
Discrepancy Project, An Update Report - Decanber 1996", available from EPRI, 

30 Palo Alto, CA. Studies by Lang and M.A. Bushey, analyzing 14 power plant test 
results, indicate that errors in emission rates using the EPA F Factor method 
could range from -8 percent to +4 percent; and, when studying five dozen coal 
samples, they found that for 18 percent of the samples the error exceeded 5 
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percent. See F.D. Lang and M.A. Bushey, "The Role of Valid Emission Rate 
Methods in Enforcement of the Clean Air Act**, EPRI 1994 Heat Rate 
Improvement Conference, May 3-S, Baltimore, MD. 

Another related art to the present invention was developed by Roughton 
5 in 1980; see J.E Roughton, "A Proposed On-Line Efficiency Method for 
Pulverized-Coal-Fired Boilers", Journal of the Institute of Energy, VoL20, March 
1980, pages 20-24. His work served in part as the basis for the above cited EPA 
methods, and is related to the Output/Loss Method. Roughton's method produces 
boiler efficiency (ilboiicr) independent of fuel or effluent flows. His work computes 

10 boil« efficiency from the process' stack and non*stack thennal losses, evaluated 
per unity of As-Fired fiiel flow. Of these losses, the major loss is the stack loss. 
Whereas this is directly measured for the Ou^ut/Loss Method, Roughton relies 
on the statistical relationship between dry effluent flow and total As-Fired fiiel 
energy flow assuming a water- free (dried) fuel. He relies on an observed ratio of 

15 these two quantities being essentially constant at 0.0008257 lbm^,,gn/BtUA5.F««iFuci 
(referred to below by the term and called the **fiiel fiactor"). The method's 
accuracy is based solely on this value remaining constant. The EPA's F-Factor 
approach fimdamentally relies on this same ratio remaining constant, see the work 
by Lang and Bushey. It has been found that for a specific fuel, having a catain 

20 Rank from a certain geographical region, this ratio is indeed constant; but found 
not the same for different fuels. Further, it makes no claim as to predicting 
heating value; indeed Roughton states: "Using this ratio it is possible to obtain the 
dry stack loss and moisture loss without the need for an ultimate analysis of the 
coal or for an accurate detennination of calorific value" of the As-Fired fuel. 

25 Cony)lete theimodynamic understanding of fossil-fired systems, for the 

purposes of improving system efficiency and regulatoiy reporting, requires the 
detennination of fiiel mass flow rate, fliel chemisto^, fuel heating value, total 
effluent flow rate, emission rates of the common pollutants, and thermal efficiency 
of die overall combustion process. All such quantities must be determined with 

30 thermodynamic consistency. There is a need for an inqnroved approach to these 
determinations. 
/// 
/// 



I 



CA 02325929 2000-09-25 
W099/4M99 PCTAJS99/0426S 

-16- 

SUMMARY OF THE INVENTION 

The iqjproach of the present invention, termed the Input/Loss Method, 
consistently determines fuel flows, effluent flows, emission rates, fuel chemistiy, 
fuel heating value and thermal efficiency, resulting in improved determinations 
5 of the thermal eflSciency of any fossil fueled system. The Input/Loss Method has 
been ^lied tiirough computer software, installable on a personal compute, and 
demonstrated being fully operational. This computer and software is termed a 
Calculational Engine, receiving data from a fossil fueled system's data acquisition 
devices. The Calculational Engine operates continuously, i.e., in "real time" or 

10 "on-line", as long as the fossil fueled system is receiving and burning fuel. 

Prior to on-line operation, the Input/Loss Method requires certain 
initiali2dng data mvolving reference fuel chemistzy and heating value, and 
reference fliel stoichiometric data associated with the reference fuel chmiistiy. 
In addition, those computer programs which will desoibe the steam generator's 

15 air handling equipment and its heat exchangers, or the gas turbine and its heat 
exchangers and/or the steam turbine cycle, require routine mitiating input. If 
operating for the first time, initial estimates of fuel chemistry and heating value 
are required, typically the reference values are used. Effluent measurements are 
required. If off-line, the assumed effluent measuremrats may be either consistent 

20 with the assumed fuel, or take a bias for the stucfy of instrumentation error. Using 
these data, oror ansiyses are preformed indicating which assumptions will 3deld 
minimum errors by exercising the Input/Loss Method as an analytical tool. Two 
such assumptions are required The first of these is how the molar fiaction of fuel 
ash, using a Moisture-Ash-Free (MAF) base, aMAF-Ath> should be treated: using a 

25 constant value of (tMAF^Aoi or correlating a^,^ as a function of MAF heating 
value, HHVm^p; measuring it directly using a fuel ash meter, or, preferably, 
determining a^,AF.Ash by explicit solution requiring the measurement of the 
system's wet combustion Air/Fuel mass ratio (a relative measurement routinely 
made in all fossil fueled systems). In general, this Air/Fuel ratio requires 

30 normalization given measurement bias such that stoichiometric consistency is 
acMeved. An alternative to the Air/Fuel ratio is the indicated fuel mass flow rate, 
also routinely measured in all fossil fueled systrais, but requiring normalization 
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(using eiTor analysis procedures of this i^^ Tte second assumption is how 

the molar fraction of fiiel wat^, using a MAF base, Of^^^AF-^ should be influenced 
by effluent water: using a constant value of the efflumt H^O, is preferred; or the 
effluent stream may be instrumented for the direct measurement of H^O. 
S When operating in real time and using the initialized data, the Input/Loss 

Method p^orms the following sequential steps: 1) obtain measurements of the 
concentrations of the conunon pollutants to accuracies conunon to the electric 
power industiy; 2) obtain measurements of the gross shaft electrical power, 
mechanical power, and/or the useful energy flow developed from the system; 3) 

10 if multiple fossil fuels are used, their properties are combined (e.g., using the 
FUEL.EXE program) to form a composite fuel, composite higher heating value, 
and, even if a single fuel, to prepare input for step 4; 4) fuel concentrations and 
heating value are input to a steam generator or gas turbine computer simulator 
(e.g., the EX-FOSS.EXE program); 5) obtain measurraients of the effluents O^, 

1 S CO2, SO2, and H2O if q)propriate, and the indicated Air/Fuel ratio (or indicated 
fuel mass flow) if q)proi»iate, to accuracies common to the electric power 
industiy; 6) the steam generator, or gas turbine, computer sunulator is executed 
producing consistent stoichiometrics given die supplied fuel and input of the 
measured effluent O2 and common pollutants, includuig the computed effluents 

20 CO2, SO2 and H^O values, the Air/Fuel ratio, the moles of fiiel per 100 moles of 
diy gaseous effluent (termed x), and at least the following thermal performance 
parameters: fiiel and efflumt flows, system thermal effldracies and emission rates 
- all consistent with the input fuel's chemistry and heating value; 7) by solution 
the molar MAF fractions of fuel carbon, water and sulfur are computed as explicit 

25 stoidiiometric solutions, both dependent principally on the effluents Oj, COj, SOj 
and H2O (which was not done in the *470 patent or the *420 patent); 8) through 
depCTdency on the molar MAF fraction of friel carbon, the molar MAF fractions 
of fuel nitrogen, o^^gen and hydrogen are determined; 9) as optioned in the 
initialization, the molar MAF fraction of friel ash is determined, the preference is 

30 by e^licit solution; 10) the molar MAF fuel species are converted to a molar dry 
base, then converted to a molar wet base (As-Fired), and finally to wet weight 
fractions (As-Fired); 1 1) the hi^er heating value is computed based on changes 
to the fuers MAF constituents, then convoted to a dry base, and then to an As- 
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Firedbase; 12) the results ofthe last two stqis, fuel chemisbya^ 
are then input to the FUEL£XE program (or a similar program) of step 3 and the 
processes repeated until convergences on the fuel moles (x), As-Fired heating 
value (HHVaf). fuel flow and minor stoidiiometric terms are achieved; 13) after 
5 convergOTces are achieved, and if errors are within criteria, the procedure is 
deemed successful, if not instrumentation is identified by the process allowing for 
both correction and the minimization of errors through application of 
multidimensional optimization tedmiques; and finally, 14) after convergences are 
achieved and error analysis completed, the operation of the system is adjusted to 

10 improve its thermal efficiency and/or to minimize the polluting emissions and/or 
to report effluent flow and emission rates to regulatory authorities. When this 
process is completed all objectives of this invention will have been met. 

The fossil flieled system operator has assurance of complete 
thermodynamic understanding of the system because of: 1) e?q)licit relationships 

1 S between measured effluents and the key fiiel constituents of carbon, water aiid 
ash; 2) an explicit relationship between these and the computed heating value; 
3) an explicit relationship between the fiiel energy flow (heating value and flow) 
and the useful »ergy flow developed from Ae combustion gases; and 4) an 
mplicit relationship between fiiel flow and effluent flow. 

20 The q>paratus necessaiy for practicing the preset invention includes any 

measurement device (or combination of devices) which determines the effluent 
concentrations of and COj,, and, if s^ipropriate, effluent HjO, CO & SOj to 
current standards found in the electric power industiy. Further, the system's 
routinely measured Air/Fuel ratio (or fuel mass flow) is required provided it is 

25 consistent, not necessarily accurate (normalization of the signal is provided) if fuel 
ash is to be determined in real time. Further, routine and conmion thermal syst^ 
data, all of which are typically known to thennal system operators, is required 
such as: effluent gaseous temperature; combustion air psychrometrics; working 
fluid flows, pressures and temperatures at key heat exchangers, and the like. 

30 It is therefore an important object of the present invention to provide a 

procedure for determining the energy flow of tiie input fuel to a fossil fueled 
system without direct measurement of the fiiel flow rate or heating value or fuel 
chenustiy, and in accomplidiing this, to assure that system efficiency, and system 
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mass and energy in-flows and out-flows are consistent 

It is a further object of the present invention to provide a quantitative 
procedure of demonstrating the consistency of the Method's results. 

It is a further object of the present invention to provide a means for 
5 determining the energy flow of the input fuel of a fossil fueled system by 
predicting the composition of the input fuel, including its ash content, and with 
this information predict heating value, and th«i back-calculate the hq)ut fuel flow 
rate from a classical use of system efficiency. 

It is a further object of the present invention to provide a means for 
1 0 determining both the total efflumt flow rate (cubic feet/hour or poimds/hour), the 
emission rates (pounds/million-BtUfue,) and flow rates (pounds/hour) of all 
effluents including the common pollutants produced firom a fossil fueled system 
by determining the fuel flow rate indirectly and having knowledge of the fueVs 
chemistry and effluent concmtrations. 
IS It is a further object of the present invention to provide a procedure for 

determining the thermal efficiency of a fossil fueled system without <tirectly 
measuring the input fuel flow rate. 

It is a fixrther object of the present invention to provide an intrinsic self- 
chedcing procedure of the Method of this invention, computed in real time, in- 
20 process, in which the computed and measured effluent COj, efflu^it HjO, 
Air/Fuel ratio and computed ftiel £actor are compared for consistent, and if errors 
are within criteria the procedure was deraied successful, or not 

It is a further object of the present invention to provide a procedure to 
identify which of the effluent CO2, efiQuent H^O and Air/Fuel ratio measurements 
25 is producing erroneous measurements, such that corrective actions may be taken 
if appropriate. 

It is a further object of the present invention to provide a procedure in 
which the errcMs made in predicting the concentration of fuel wato- will be oflf-set 
by fuel ash, though an explicit computation of fuel ash. TTius any error in fuel 
30 water, as off-set by fuel ash, will have negligible effect on the As-Fired heating 
value. 

It is a further object of tfie present invention to provide a means for 
determining the flow of a fiiePs solid non-combustible mineral material 
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(commonly referred to as fuel ash) associated with coal fuel, and thus adjust and 
improve the operations of the ash removal equipment in a combustion system^s 
effluent stream. 

It is a further object of the present invaition to demonstrate that all 
S aformirationed objectives associated with a fossil fueled system are also, and 
herein declared, objectives associated with any fossil fueled system producing 
electrical power, mechanical power and/or useful energy flow from the system. 

It is a further object to provide an approach which yields improved results 
as compared with EPA Methods 1, 2, 17, and/or 19. 
10 It is a further object to provide an approach which yields improved results 

as compared with the methods of the ASME Power Test Codes 4. 1 , 4.4, and/or 22. 

Other objects aiid advantages of the present invention will become 
apparent when the Method and apparatus of the present invention are considered 
in conjunction with the accompanying drawings. 

15 BRIEF DESCRIF nON OF THE DRAWINGS 

FIG« 1 is a block diagram illustrating the initialization procedures 
including the study of instrumentation errors; 

FIG. 2A and 2B represmt a block diagram of the present approach, 
showng the determination of fuel flow, system efficiency, fuel chemistry, fiiel 
20 heating value, effluent flow and emission rates herein collectively referred to as 
Fia2. 

FIG. 3 is a plot of Moisture-AslinFree molar fuel hydrogm versus 
Moisture-Ash-Free molar fiiel carbon using typical coal data 

FIG. 4A and 4B represent a table of suggested Calculational options 
25 associated with this invention allowing for the analysis of any hydrocarbon fuel 
and coals of various Ranks herein collectively referred to as FIG.4. 

FIG. 5 is an mor plot showing laboratory-reported wet heating values 
versus those computed by the Method of this inv^ion for a coal-fued power 
plant having wildly varying fuel chemistries and heating values. 
30 FIG. 6 is a schematic representation of a conventional or fluidized bed 

power plant illustrating use of stoichiometric relationships. 
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DESCRIPTION OF THE PREFERRED EMBODIMENT 

For clarity, note that FIG. 6 depicts fossil-fired pow^ plant thermal system 
20. In this power plant thermal system 20, a fuel feed 22 and combustion air 24 
are provided finom an upstream side 26 of the heat exchangers/combustion region 
S 28. Note that this region 28 does not include the air pre-heater 36. In addition^ 
in some types of power plant th^mal systms 20, other materials may be ii^ ected 
into the heat exchangers/combustion region 28, such as a flow of limestone 3 1 to 
minimize effluent SO2 by chemically binding sulfur. The fuel feed 22 contains, 
in general, combustible material, water and mineral matter (called fuel ash). The 

10 &el ash is an unbumable component that passes through the system with little 
physical change, but which is heated and cooled. 

In the heat exchangers/combustion region 28, the fossil fuel 22 is burned 
with the combustion air 24 to form hot combustion products. Heat from the 
combustion products is transferred to a working fluid 30 that flows through heat 

IS exchangers 32 that are depicted as integral with the heat exchanger/omtbustion 
region 28. The heated workiqg fluid 30a is used in a manner appropriate to a 
working fluid to generate a useful ou^ut 33 (for a conventional power plant such 
usefiil output may be supplied to a turbine cycle which could produce electricity). 
There may be a water leakage 29 into the hot combustion products, not associated 

20 with the fiiel feed 22 as a result ot for example, soot blowing or tube leaks fix>m 
the heat exxAmgers, 

After leaving the heat exchangers/combustion region 28 on its downstream 
side 34, the cooler combustion products commonly flow dirough ducts 35 which 
may contain effluent ash removal equipment, passing then to an air pre-heater 36, 

25 where a further portion of their heat energy is transferred to an incoming air 
stream 38, v^ch Aea becomes the combustion air 24. The total air delivered to 
20 is 25. In many cases, a leakage flow 40 of air enters the flow of combustion 
products from the incoming air flow 25 as it passes through the air pre-heater 36. 
The fiirthCT-cooled combustion products leave the air pre-heater 36 and pass to a 

30 stack 42 and are then exhausted to the aUnosphere. 

This gCToal system description is applicable to a wide variety of fossil- 
fued power plants, such coal-buming power plants, oil-burning power plants, gas- 
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fired pow^ plants, biomass combiistors, fluidized bed combustors, a conventional 
electric power plant, a steam graerator, a package boiler, a c<xnbustion tuibine, 
and a combiistion turbine with a heat recoveiy boiler combustion tuit This 
list is not meant to be exhaustive, however, and is presented to illustrate some of 
5 the areas of applicability of the present invention. 

It is desirable to quantify the operation of this powa- plant thermal system 
20. If it is characterized quantitatively, then the relationship of variations in the 
operating parameters to the inpxsts and outputs of the power plant thermal system 
20 may be understood This understanding, in turn, permits the operation of the 

10 power plant thermal system 20 to be optimized for thermal efficiency and 
pollution minimization. 

In this system, some quantities are readily measured, and others cannot be 
measured on-line, in real time with accuracy sufficient to quantify the operation 
of tiie power plant thermal system 20 to the required accuracy. For example, gas 

1 5 concCTtrations, working fluid flow rates, and temperatures and pressures may be 
readily measured with good accuracy convmdonal sensors located at a defined 
systrni boundary 44 and elsewhere sudi as in the regions 35 and 42. Ontheother 
hand, the As-Fired fuel chemistry of the fuel feed 22 cannot be measured 
accurately m real time if the fuel is coal or other bulk fuel. The present invention 

20 uses a modeling analysis to quantify die As-Fired fuel chemistry of die fuel feed 
22 and other param^ers in a self-consistent manner, which allows the entire 
power plant thmnal system 20 to be quantitatively understood. Once the power 
plant syston is quantitatively und^-stood, its operating param^ers may be varied 
to optimize its thermal performance. 

25 The present invmtion provides an Input/Loss Method which allows for a 

complete understanding of fossil fueled combustion systCTis such as power plants, 
through application of non-direct but explicit determination of fuel and effluent 
flows, fuel chemistry, fuel heating value and thermal efficiency, resulting in 
improved tfiermal efficiency. In addition the Input/Loss Method of the invention 

30 determines the effects of inaccuracies in the Method's primary measuremaits. 
This invention was developed primarily for coal-fired systems. However, this 
invention may also be implied to any fossil fueled combustion process including, 
but not limited to, gas, oil, bio-mass, slurry fiiels, wood, agricultural byproducts 
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sudi as shells from nuts, trash and refuse. 

More particularly, this invention relates to the Iiq)ut/Loss Method for 
monitoring the operation of a fossil fueled process by analyzing the composition 
of combustion eflQuents CO2, 02, SO2 and, deprading on the characteristics of the 
5 fuel and combustion process, effluent HjO and/or the systems' indicated 
combusticm Air to Fuel ratio. Other routine ^stem data are also required such as 
gaseous effluent temperature. Having computed the thermal efficiency and 
consistent fuel flow rate, and having computed the fuel's chemical composition, 
and having computed the fuel's heating value, the plant's effluent flow rate and 
10 emission rates may then be determined. Further, this invention allows for the 
determination of tfie effects of inconsistencies in the principal measurements of 
effluent COj, effluent HjO and the systems' Air/Fuel ratio (or indicated fuel flow), 
as sudi inconsistencies impact the thermal efficiency of the system and its effluent 
flow rate. 

1 5 The subject of this invention includes any fossil-fiieled system producing 

a useful and measurable output Such fossil-fossil systems include, for example, 
conventional electric power plants; fluidized bed electric power plants; any boiler 
whose combustion gases cause tiie increase in energy of a measurable fluid; a 
combusticm turi>ine producing electrical or mechanical power; and a combustion 

20 turbine producing both electrical and/or medianical pow^, and whos^ combustion 
gases cause the increase in energy of a measurable fhiid. Although the traditional 
working fluid from fossil-fired systems is steam, this invention places no 
limitations on the nature or quantity of useful output associated with fossil 
combustion. However, if fuel and effluent flows are to be determmed, the 

25 invention does require that the working fluid's energy received from the 
combustion process be measurable. This invention applies to a wide range of 
systems, from small boilers as might be used for space-heating hotels, to large 
commercial electric power plants. Working fluid flows, such as steam flow, could 
typically range from 10,000 pounds/hour (Ibm/hr) to 10,000,000 Ibm/hr. 

30 Electrical or mechanical power could typically range from 5 Mega- Watts (MW) 
to 1,500 MW. 
/// 
/// 



I. 
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Svstem Calculations 

The present invention allows for a complete understanding of fossil-fueled 
combustion systems through explicit determinations of fuel and effluent flows, 
fiiel chemistiy, fuel heating value and ttiermal efficimcy fluou^ thermodynamics, 
5 not through direct measurement of fuel flow, fiiel chemistiy nor through effluent 
flow. The approach principally relics on measuremcaits of the efilurats from 
fossil-fiieled systems, and otho- routine syston data. Given the nature and choice 
of such data, it has the potential of high resolution on a contmuous basis, but the 
process only anticipates common industrial accuracy standards of required 

1 0 instrumentation. The data may be input to a computer program for resolution of 
mass and energy balances associated with the system. Measurements include the 
concentration of combustion gases exiting the system (i.e., smoke stack), the 
direct production of power, and the total energy deposition to the working fluid 
from the hot combustion gases. 

1 5 The thermal efficienQr of a fossil-fueled system, termed ^stem diermal 

tfRdency is defmed in the usual manner (Le*, total useful output divided by total 
inputs): 

n.^ — Gross Shaft Power fCn T Jsefiil Hnergv Flow Delivered 

Fuel Energy Flow -f System Energy Flow Credits (1) 

20 In Equation (I), "Gross Shafl Power (CT)" refers to direct shaft power developed 
from the combustion gases as normally associated with a combustion turbine. 
Such a system may or may not CTiploy a heat recovery steam generator (if it does 
then the "Useful Energy Flow Delivered" term will be non-zaro). For a 
conventional fossil- fueled power plant the Gross Shaft Power is zero, and the 

25 Usefiil Energy Flow Delivered refers to the raergy flow to the working fluid of 
a turbine cycle (TC). In general the terms Gross Shaft Powct (CT) plus Usefiil 
Energy Flow Delivered define the thermal system's "useful output". For a 
convoitional power plant or a conventional boiler, it is necessaiy to separate the 
powffl- production process, or turbine QTcle, from the combustion process. Note 

30 that the efiFiciency of the power production process could be thought of as the 
efficiency associated with the working fluid's process. 
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^turtw^cyde ^ Useful Powcf Produced CrC\ 

Useful &iergy Flow Delivered (2) 

Therefore: 

'Isystcra ^ ^boiler ^tiiitNiie<yde (3) 

where for a conventional boiler: 



^boacf= : — Usgfal Engrey Flow Pelivgrgd 

Fuel Energy Flow + System Energy Flow Credits (4) 

For a combustion turbine system Ti^rtm«ydc is unity, and the term is 
identically ti^^^ as defmed by Eq. (1). In all further discussions the term Tiboi,^ 

10 is used in place of t]^^ v^dien referring to a combustion turbine system 
(understanding that the value of ti^^^^^s loO), and the "Useful Energy Flow 
Delivered" of Eq.(4) is replaced with "Gross Shaft Power (CT) + Useful Eneigy 
Flow Delivered" of Eq. (1). Further, for convenience of discussion, combustion 
tmbine efficiencies employ a lugger heating value; conversion to a lower heating 

1 S value base, as commonly used for combustion tutbine systmis in North America, 
is a simple conversion ttiBt is well known in the art and described in ASME PTC 
22. Outside of North America, lower heating values are used for most analyses 
of fossil fueled systems. Several authors in this field disagree with the use of the 
tenn "System Energy Flow Credits"; see the denominator of the efficiency 

20 EquaticHi (4). Its use, or not, is not important to this invmtion; however, its 
consistent treatment relative to tiboua is most important If energy credits are to be 
used, which is the Prefiored Embodiment, then the tmn must appear in the boiler 
efficiency calculation, and also when fuel flow is computed, i.e., Eqs.(19B) and 
(21). If the term is not used, thm only fuel flow and heating value will define fiiel 

25 enei^gy flow and thus T]boiicr • 

The evaluation of electrical or mechanical power terms is well known. The 
Preferred Embodiment uses gross electrical or mechanical power, that is total 
power derived directiy from the combustion gases. However, net power may be 
used, that is gross power less that power consumed within the system, called 

30 ^'house loatf\ Use of either of tfiese terms is not important to this invwtion, only 
that gross or net terms be used consistently. The evaluation of the term Useful 
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Energy Flow Delivered is well known, and involves simply measuring of mass 
flow rate, and pressures and tonperatures of the working fluid. Therefore the 
evaluation of ntt^rt,ine<ydc is known to those skilled in the art, and may be 
accomplished with high accuracy and reasonable expense, 
5 It is evident fix)m the above discussion that all principal terms are 

descriptive of any fossil*fueled system producmg useful energy flow and/or power 
(in any form). The Mettiod of this inv^tion is equally applicable to a combustion 
turbine as to a conventional boiler producing useful energy flow by heating a 
working fluid such as steam or liquid water, and to a conventional fossil-fueled 

10 power plant consisting of a conventional boiler heating liquid water to create 
steam for direct usefiilness and/or for a steam cycle producing any combination 
of electrical or mechanical power, defined as useful power fiom a turbine cycle. 
The Method of this invention is equally applicable when usmg higjier or lower 
heating value as a bases for the efflciency definition.. 

15 The boiler efiiciency of Eq. (4) is helpful for continuous monitoring of 

thermal performance if the fuel flow and other parameters comprising may 
be measured accurately. Even if so monitoring, then the directly measured fuel 
flow, fuel wat^, fuel ash and/or the resultant boiler efficiency may be compared 
to the calculated values via procedures herem disclosed. The Method of this 

20 invention computes altfiough Eq.(4) (or Eq.(l) for a combustion turbine) is 
reformulated sudi that fiiel flow is caiculationally excluded. Aft^ excluding the 
fuel flow term, three majcM- deficimcies in the knowledge of a fossil-fiieled 
system's thermodynamic process remain, and especially vAien considering a coal- 
fired system: 1) the complexities of the combustion process (and if a combustion 

25 turbine, the complexities of the compression and expansion of gases); 2) the 
specification of thermal losses not directly related to the combustion process (but 
which could affect measured fiiel flow); 3) the variability or uncertamty of the 
heating value; and 4) the variability or uncertainty of the fuel chemistry as it 
impacts Tiboijg^ and the heating value. 

30 A computer program EX-FOSS.EXE has htcn developed to address the 

first two of these difficulties. It is a commercially available program which has 
been in use in ttie power generation industry since 1985 and is available (together 
with the I^L.EXE program) fix)m Exergetic Systems, Inc. of San Rafael, 
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Califomia. The methodology of the EX-FOSS.EXE program is presented in detail 
in the *420 patent, see its Figures 1 A and 1 B, and their associated discussion. The 
FUEL.EXE program, using both known fuels (e.g., stabilizing gas, oil and/or 
coals) and the fiiel being resolved, performs simple weighting calculations 
S resulting in a composite fuel. FUEL.EXE may accept any number of fuels, each 
of which is described by . its elementary molar or weight fractions. Gaseous 
components are input by molecular molar or weight fractions. Each descriptive 
filers fractions must sum to unity. These are combined to produce a composite 
fuel, at the user's option, using a molar or weight base. The individual fuel's 

10 heating values are also weighted. This composite is then supplied to EX- 
FOSS.EXE as required input. Hie EX-FOSS methodology separates the 
definition of boiler efficiency into components which, taken separately, 
calculationally exclude the first of the aforementioned problem areas. When 
separated, terms called combustion efficiency, r\c^ boiler absorption 

15 efficiency, Tj^, are developed. Use of rj^ allows consideration of only input and 
output terms which may be measured with high accuracy, thus eliminating 
resolution of combustion's thermodynamic complexities. t|c addresses the major 
loss in a fossil fired system, i.e., stack losses. The use of allows for the 
consideration of only "non-stadc" losses, losses viiich generally have minor 

20 impact on system efficiency, thus whose understanding does not require high 
resolution. The following terms are defined, and typical units are provided for 
clarity: 

m^F s As-Fired Fuel Mass Flow Rate; IbmAp/hr. 
HHVP H Higher Heating Value at Constant Pressure; Btu/lbm^f. 
25 HBC 5 Specific System Energy Credits; Btu/lbm^. 

EF s Fuel Energy Flow (fuel flow x higher heating value 
at constant pressure); Btu/hr. 
sm^^pHHVP 

CF = System Energy Flow Credits (fuel flow x specific energy 
30 crecUts relative to the system); Btu/hr. 

s m^^pHBC 

HPR s Enthalpy of the Combustion Products (includes the heat of 

formation plus /CpdT at the effluent boundary); Btu/lbm^^p. 
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HRX s EtHbaply of the Reactants (based on die heating value, 

sensible heating and eneigy credits); Btu/lbmAp. 
HSL s Specific Stack Losses (includes losses directly 

affecting the oiergy released during combustion, defined by 
5 PTC4J:Lo,L„^L„,L^,Lx,Lz,Lco,Lu„&Luw:;all 

divided by m^); Btu/lb>vF. 
HNSL = Specific Non-Stack Losses (boiler or combustion turbine losses 
whose medianisms originate fi^m the combustion process 
or fi"om the hot gases, and interface directly with the 
10 environment thus a direct effect on the BBTCtem; defined 

by PTC 4.1 as: Lj), Lp, La, L, & L^c; defined by PTC 4.4 as: 
Lfl & Lw; in addition, if appropriate, to the gas turbine's 
shaft bearing losses, applicable gear losses and other similar 
finctional losses associated with deliveiy of useful shaft 
1 5 power directly by combustion gases (all divided by mj^); 

Btu/lbmAF. 

ERC = Enetgy Released during Combustion; Btu/hr. 
= EF + CF - E(Stack Losses) 
= EF + CF-mA,HSL 
20 =mA,<HPR-HRX) 

BBTC = Gross Shaft Power (CT, if a combustira turbine) 
plus Useful Energy Flow Delivered; Btu/hr 
= ERC - £(Non-Stack Losses) 
= ERC-mAfHNSL 

25 With these definitions, equivalent ways to esxptess fossil boiler efficiency include 
the following. 

Ibone,-- BETCL 

EF + CF (5) 



30 



n^n— ERC - ErNon-Stflcfc I^ses;> 
HF + CF 



(6) 
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tlK..^, fEF + CF^ - EfStack Tx>sses^ - EfNon-Steclc T^^q^;) 

EF + CF (7) 

EF + CFl ERC / (8) 

S £q.(8) suggests that boiler efficimcy may be divided into two separate 
efficiencies: one descriptive of the combustion process per se (called the 
combustion efficiencyX and the other descriptive of certain non-stack losses 
(called the boiler ahsoiption efficiency). 

The combustion efficiency definition is suggested from the Input-Ou^ut 
1 0 Method defined in PTC 4. 1 or in PTC 4.4: that is, net weigy flow produced at the 
thermoc^iiamic boundary divided by the total enetgy flow input (fuel en^:gy flow 
and system eneigy flow credits). The following develops rfc on a unity fuel flow 
bases: 

rie^ 1.0- !S(Stapk Losses) 
15 EF + CF 

= (EF -^ CF) ' gfStack Imm) 

EF + CF (10) 



EF + C 

20 Tic = [mAp(HPR-HRX)]/ (EF + CF) (12) 



EF + CF (1,) 



= HPR* HRX ■ 

HHVP + HBC (13) 

The boiler absorption efSciency is derived from the Heat-Loss Methods 
found in the PTC 4.1 and 4.4 although limited to non-stadc enei^ temis. It miist 
25 be ref4»imced to the Bietgy Released during Combustion term (ERC) if all losses 
(on a systems basis) are to be additive when calculating the boil^ efficiency, see 
Eq.(7): 

1,0, £(Non-Stflck Losses) 

T]c(EF + CF) (14) 
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The quantity iic(EF + CF) defines the ERC tenn, see definitions and Eq.(l 1), 
thus: 

iIa=1.0 - niAiHNSL/ERC (15A) 
Ti^==1.0- HNSL 

5 HPR-HRX (I5B) 

The Stack Losses includes the following PTC 4.1 terms relating Stadc 
Losses to total As-Fired fuel flow rate: 

mAFHSL = Lo-i-L^ + LH + L„^ + Lx + Lz + Lco-»-Lo„-hLu„j. (16) 

The quantity HNSL includes die following PTC 4.1 tenns, PTC 4.4 terms, 
10 and terms relating to turbine losses (energy flow losses in delivoing gross shaft 
power): 

m^iHNSL '=Lft4.i + li,+Lu+Lr + lw + U4.«+lw + (turbine shaft losses) 

(17) 

The combination of the combustion efficiency and boiler absorption 
1 5 efificioicy is equivalent to the classical definition of boiler efficiency as defined 
by ASME Power Test Codes (e.g., PTC 4.1), British Standards (e.g., BS2885), 
and othos. The following, using direct energy flow terms, again demonstrates the 
derivatimi of boilo* efficiency, see Eq.(13) and £q.(lS): 

nboikr=ncnA (i8) 

20 Ti,,ii«,= HPR-HRX . [HPR-HRX-HNSL]/(HPR-HRX) (19A) 

HHVP + HBC 

^twilet " -^^^ " HRX ~ HNSL 

HHVP + HBC (192^ 
^ 

25 m^iOfflVP-i-HBC) (20) 

Eq. (20) may be reanan^ to solve for fiiel flow rate, provided heating vahie and 
BBTC are accurately known. The resolution of fuel flow and heating values are 
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objectives of this invention: 

in^= EBI£ 



Tl^a,XHHVP + HBC) ^21) 

By separating boiler efBciency into combustion and boilo* absorption 
S components, the analyst has knowledge as to where degradations are occurring. 
If combustion efficiency decreases (i.e., stack losses increase), the plant engineer 
would consider fuel-air mixing equipment, degradation of hardware directly 
interfaced with delivery of the combustion air, low heat content in the fuel, 
improper operaticm of the fuel's burner mechanism, improper operation of the 

1 0 combustion proper such as tiie location of the fu^-ball, etc, - all sources directly 
affecting the combustion process (i.e., stack losses). The terms comprising 
combustion efficiency may be easily reduced to a unit basis of total as-fired fuel, 
refer to Eq.(l 3); as such these terms have the potential to be determined with great 
accuracy. HHVP is the higher heating value; corrected for an assumed constant 

IS pressure combustion process, if necessary, using the procedures established in 
PTC 4.1. 

In a similar manner, if the boiler absorption efficiency decreases (non-stack 
losses increase), consideration should be given to terms affectmg this efficiency: 
radiation & convection losses, gas turbine bearing losses if extracting power 

20 directly from tfie combustion gases, heat exchanger water/steam leaks, heat 
exchariger effectivmess, etc. The boiler absorption efficiency also has the 
potential to be determined witii high accuracy. Generally this term qsproadies 
unitj^ thus its mot is no greala- than its compliment (if = 98 percent, one miay 
safely assume a maximum error of ±2 percent). Altiiough is dq)endent 

25 (through flie torn ERC) on ifc; and a given degradation m rjc will effect tIa, the 
impact on relative changes is general!/ small Also, by an iteration technique, T)a 
may be resolved without a priori knowledge of fijel flow rate. Thus, both Tfc and 
tJa, therefore r\y^^ may be determined indep«ident of fuel flow. 

The enthalpy of the products (HPR) may be accurately calculated using 

30 methods discussed in the '470 patent. The enthalpy of the reactants (HRX) may 
be accurately calculated using metiiods discussed in die '470 patait, provided the 
eneigy associated with injected limestone is included, if applicable, using sunilar 
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techniques. 

The Plreferted Embodimem employs water and st^^ 
internally consistent within ±0.0008% maximum for enthalpy calculations. They 
derive from the 1969 woric by Keenan, Keyes, Hill and Moore, modified by 
S NASA. Combustion and atmospheric gases have total consistency (zero error) 
given their explicit solution techniques between heat capacity, enthalpy and 
entropy; they are based on American Petroleum Institute standards. Moist air 
psychometrics are consistent within ±0.0010% for common enthalpy calculations. 
Below 38F, as low humidities are approached, consistency is within ±0.0090%. 

10 Heats of Formation are obtained from the world-recognized CODATA standards 
for thermodynamics. 

The above description of system equations and related procedures results 
in the concept that ftiel flow rate may be determined with complete consistency 
from thmnal efficiency. The accuracy of the computed fuel flow is intrinsically 

15 dependent on fiiel chemistry and heating value, addressed below. The 
computation of boiler efficiency using any operational and consistent method 
available is acceptable, provided it is independent of fiiel flow such that fuel flow 
may be back-calculated from the classical definition of thermal efficiency, and 
thus is thermodynamically consistent with Usefiil Enogy Flow Delivered from the 

20 system (and/or shaft power delivered from a c(»nbustion turbine). As used by the 
Input/Loss Method, the determination of boiler efficiency requires only 
computational consistency, not absolute accuracy. Computationally, boiler 
efficiency need only produce an accurate fuel flow, consistent with energy flows 
delivered BBTC(Eq.(5)). 

25 Combustion Equation and Fuel Chemistrv 

Whatever specific procedures are en^iloyed to determine boiler efficiency 
(or system efficiency for a combustion turbine), this invention is based on a 
consistent treatmrat of the combustion equation. Any determination of efficiency 
is adequate when applied for this invention, provided the terms of the combustion 

30 equation are treated in an operational and consistent manner. The following 
combustion equation is taken from the EX-FOSS computer program for 
illustration of the preferred procedures. However, another consistent and 
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available method for tfie detenniiiation of efficiency may be found in Chapter 9 
of Steam. Its Generation and TTse^ 40th Edition, edited by S.C. Stultz and J.B. 
Kitto, published by the Babcock & Wilcox Company . Barbeiton, OR As applied 
in the present invention, ''consistmcy'' is defined by three attributes: 1) use of 
S accurate molecular weights and consistently used thermodynamic properties; 2) 
use of accurate and correct molar balances, proof of which is demonstrated when 
such molar balances are used as a bases of system mass flow balances yielding 
minimum enx>r; and 3) the incorporation of these same molar quantities within the 
computation of boiler efficiency. This later attribute is most important to the 

10 present invention. The concern is not the denominator of the boiler efficiency 
equation and whether boiler credits are included or not, but rather the equation's 
number. The terms HPR and HRX discussed above, as integrally based on the 
molar quantities defined below, with HNSL, defmed the usefiil power and/or 
energy flow delivered by the system per unit of fuel; thus use of these terms 

15 demonstrate required consistency. The equality mAi^HPR - HRX - HNSL) = 
BBTC is part of this consistency. 

In the present approachi the As-Fired fuel chemistry is first calculated 
based upon system measurmients and the use of a mathematical description or 
model of the diermal system. All calculated system parameters (e.g., thermal 

20 efficiency, total effluent flow, emission rates, fiiel flow, heating values and 
parameters associated with Second Law analysis) are thereafter calculated 
principally fiom the As-Fired fuel chemistry, and are therefore necessarily 
''consistent with each other. By contrast, in each prior approach the various 
syst^ parameters are either not calculated at all, or calculated using a variety of 

25 , calculational techniques and with diffa:ent input information, with the result that 
the system parameters typically are inconsistent with each other to some degree. 

The preferred stoichiometric equation relating reactants to products is 
presented as Eq.(29). The nomaiclature used in Eq.(29) is unique in that brackets 
are used for clarity: for example, the expression "OjlHjO]'' means the fuel moles 

30 ' of water, algebraically simply a^; the expression "d[COJ" means the effluent 
moles of CO2, algebraically simply "d**. Note that molecular quantities are used 
versus atomic for nitrogen, oxygen, and hydrogen; conversion is routine. The 
quantities comprising the combustion equation are traditionally based on an 
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assumed 100 moles of dry gaseous product. This assumption is useful when 
measuringstack«nissionssincetheconunonlymeasuredconcenlrad(Hisarebased 
on diy molar fractions; if measured on a wet base, such measurements may be 
converted to a diy base. 

5 xIooCCyrHzJ + a,[NJ + OitHjO] + a,[OJ + 04[C] + ajtHJ + aJS] 

+ o,[COJ + o,[CO] + o,[H2S] + a,o[ash]}^^p^ 
+ bJHPl,,^ + {(1 + p)(a[OJ + a(|>^[NJ + b^[HjO])U 

PLS 

= d[COJ + g[OJ + h[NJ + j[HjO] + k(SOJ 
10 + {e[CO] + £[HJ + /[SO3] + m[NO] + ppSTjO] 

+ qlNOj] + t[CYPiHzp,] + u[CYp2H2P2]}Minor Components 

+ obpi^[CaS04-zH20] + {(1 - 0 + Y)bpLs[CaO]}E«essPLs 
+ xa,o[ash] + v/[Cj^^ + "^[C^r^] 

+ {pafOJ + MaJP^i] + PbAEHjOlU^ (29) 

1 S The following is presented to assist in understanding Eq.(29), where all are molar 
quantities unless otherwise indicated: 

X =Moles of As-Fired fuel per 100 moles of dry gas product (the solution base), 
aj - Fuel constituents per mole of As-Fired fiiel (wet-base), Ea^ = 1 .0. 

ttj are derived from weight fractions as input to FUEL.EXE. 
20 hj^ == Moles of moisture in entering combustion air, per base. 

b^ = 0)(1.0 + 4>Act)^Ai/NH2o» where o) is specific humidity of the local 
atmosphere determined from psychrometric measurements. . 
pby^ = Moles of moisture entering witii air leakage, per base, 
bz = Moles of water/steam in-leakage from the working fluid, per base. 

25 b^ =xN^E[m,«,^(mAFNH2o)l 

bpLs = Moles of Eure iim^tone [CaCOj] required for zero 

CaO production, per base. 
bpLs = - k^a ; where: kp s x(a6 - a^) / [1.0 + {l/k^),^, v/here the 

ratio of SO3 to fuel generated SO2 (i.e., before reaction with 
30 limestone) (t/kf)iapuiy supplied input, is the measured 

effluent SO2, without limestone injected, kp = kj^. 
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Molar ratio of excess CaCOj to stoichiometric 

CaCOj (e.g., Y = 0.0 if no effluent CaO). 
(niuj/mAF)x Naf / (5 bpuj Ncacw) - 1-0; where m.^ is the system's 

indicated actual limestone mass flow rate, and ^ is a 

mass ratio of actual limestone to pure CaCOj it contains. 
Moles of HjO per CaS04, stpplied as input based on 

periodic laboratory analysis of boiler refuse. 
Kronecker function: unity if (Ofi + a,) > 0.0, zero if no sulfiir 

is present in the fuel. 
Air pm-heater dilution factor, a ratio of air leakage to true 

combustion air, a molar ratio 
m{R^ - 1.0) / [a/f^(l-0 + <|)aJ]. 

Ratio of total moles of dry gas from the combustion process before 
entering the air heater, to the total moles of dry effluent (after 
leaving the air pre-heater); the air pre-heater leakage factor 
typically obtained from either the ratio of COj or Oj across the 
air pre-heater. 

Moles of combustion Oj, per base. 

Ratio of non-oxygen gases (nitrogen and argon) to oxygen in the 

combustion air local to the system. 
(1.0-^^/^^. 

Concentration of in the combustion air local to the system. 
Molar quantities of drv gas product at the system's boundaiy for 

specific compounds: d, e, f, g, h, k, /, m, p, q, t and u. 

Hie sum is denoted as £nj. Dry air gases associated with 

system leakage (Pa and pa<|>Act) ^ treated separately. 
Molar quantities of non-gas iMOduct at the system's boundary for 

specific compounds: j, obpLs> (1 - o + Y)bpLs» xa,o, w and v. 

The sum denoted as Diji. Moisture from air leakage (Pb^) 

is treated separately. 
Molecular weight of compound or element], defined by the 

1983 lUPAC to six significant figures. 
Molecular weight of the As-Fired (wet) fiiel. 
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AF Ratio of mass flows of combustion air to As-Fired fiiel, based on 

mpui 

indicated, uncaiibrated, plant data, unitless. 

V = I(»7^Botto«^Ca>*<«>^ttom-Ash) 
5 + (HT^ny^Caibon/Fly-Ash)] xHu/Nc. WFx are weight fractions. 

Resolution of Eq.(29) proceeds by solving for all n; and nji quantities. 
Minor component terms of Eq.(29) are typically resolved either through direct 
measurement (e.g., for CO), or assumed zero values, or tiirough assumed 
relati(Miships (e.g., SOj is a fimction of an input molar ratio of SO, to SO,, defmed 
10 above. All Minor Components typically have only low parts-per-million 
concentrations thus have little impact, witii tiie possible excqjtion of the terms v 
& w. The true importance and fimctionality of Eq.(29) to the Input/Loss Metiiod 
lies in the fact that consistency of molar balances is needed for successful system 
understanding as discussed above, for conservation of mass flows, and for 
15 resolution of fuel chemistry. Fcfr clari^ the following major terms are associated 
with the system's thermodynamics: 

Total effluent water s J = j + bAp 

Total efOuent oxygen s G = g + aP 

Total efllumt = Xjij + En^j + P(a + a<|>Aet + b^) 
20 Total effluents before tiie air pre-heater =RfJ^ -^R/J^ji^. 

Total As-Fired pure limestone = (1 + Y)biu 

Diy combustion air without air pre-heater leakage = (a + 

Dry air from air pre-heater leakage presort in effluent = p(a -i- a^^^ 

Wet combustion air * (a + a4)Aa + b^) 
25 Total wet combusticm air and air pre-heato* leakage 

= (l + P)(a+«l>Act+bA). 



Eq.(29) is unique in describing at least three features of critical importance 
when improving die perfonaMnce of conventional power plants witii the Metiiod 
of tills invention. These include: 1) its abUity to address air pre-heater leakage 
30 tiirough application of j?Act through tiie dilution term P; 2) tiie ability to describe 
effluent concentrations on eittier side of tiie tat pre-heatCT, again tiirough 
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^plication of R^, 3) its ability to address limestone injection commonly used in 
fluidized bed combustors; and 4) the use of an explicit term allowing for 
variable O2 concmtration in the system*s local combustion air. Air pre-heater 
leakage dilutes all combustion effluents with moist air from the local environment, 
S thus all important effluents H^O, CO^ and O2 used for this invention are altered. 
Furthermore, many times a power plant's more precise effluent measurements, 
especially O2, may be found on the air pre-heat^-'s inlet (economize outlet), and 
not at the air heater outlet; refer to FIG. 6. Although most environmental 
regulation requires effluent measurements at the system's boundary, translation 

10 between the air heater inlet and outlet measurements is many times essential. 
Eq.(29) allows for such translation through the term, defined above such that 
1 00 moles of dry gas are computed both at the upstream and downstream locations 
of the air pre-heater, see ''Boiler" of FIG. 6. Thus effluents may be used by the 
present invention either upstream or downstream of the air pre-heater; refer to the 

IS G and J terms defined above, allowing conversion between measurements with 
and without air leakage. Combustion gas conditions upstream of the air pre-heater 
and before the exit of the heat exchangers and combustion region, see FIG. 6, 
would employ terms of gR^ and ]Rj^ Combustion gases downstream of the air 
pre-heater typically exit the system to the environment, in other words the gaseous 

20 effluent boundary of the system. Limestone injected into the combustion process 
will create additional effluent CO2, and could decrease die effluent Hp if the 
sulfete product is malrixedwitti water, CaS04zH20. Of course COjand HjOare 
two important effluents to the present invention. In addition to ttie basic 
stoidiiom^cs associated with limestone injection, the detail afforded Eq.(29) 

2S allows numerous and obvious determinations of molar and mass ratios. For 
example the molar ratio of fiiel sulfur to limestone's calcium, termed die 
limestone conversion ratio, important to a fluidized bed operator, is simply k^[(\ 
+ Y)bpLsl- The concentration of Oj in the system's eombusiion air (termed A^, 
leading to the tmn 4>Act) affects all system effluent and combustion air terms. 

30 In summary, the aforementioned technique describes the process of 

calculating boiler efficiency of both any thermal system and in particular of 
practical conventional power plants based on effluent measurement data, fuel 
heating value and several parameters of minor importance. The next stage of the 
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process involves the recognition that a given fuel has an unique chemical 
composition^ thus when burned will yield unique concentrations in its gaseous 
efiOumt The principal gaseous effluents fipom any fossil combustion process are 
N2, CO2, H2O, O2 and SO2. H2O, vAim effluent from combustion, is in its super- 
S heated phase, thus acting as a gas. The source of N2 is principally the air used to 
bum the fiiel and it has littie chanical reacti veness, thus its sensitivity to the fiieFs 
chemical composition is not significant SO2 effluent concentrations are generally 
in the parts per million, thus its impact has minor importance. 

The relative concentrations of carbon, a^j and hydrogen, found in any 

10 fossil fiiel will have significant impact on the relative concentrations of CO2 and 
H2O found in the effluent. In addition, these effluents will be influenced by the 
following: Oj used to bum tiie fuel (i.e., the Air/Fuel ratio); fuel water, in- 
leakage of water; and water in the combustion air. This implies that the molar 
fractions of CO2, HjO and Oj present in the effluent (the system's boundary, i.e., 

15 its smoke stack or translated from data from the air pre-heater inlet) must be 
unique relative to the fuel input and supplied combustion air. 

The following elementary molar balances may be derived from the 
combustion Eq.(29). The expressions are simply convenient groupings of 
quantities, principally comprising measured effluents (known values) which have 

20 the greatest influence on the individual element of interest. Nfinor fuel terms 
carried within expressions are multiplied, initially, by an estimated fiiel moles, 
X. These mincM-tcrnw are qiuckly resolved when converging on X. Anexception 
is T'y^, discussed below. (Given these groupirigs the expressions of Eqs,(36) 
thm (40), with solution of the term ""a^ as discussed below, may be treated as 

25 known quantities. The elraientary wet fuel compcmmts typically associated with 
coal are considered unknowns, as are the fuel moles, these include the following: 
tth <^4^ ^5* «6» ^10 and x. Many coal-fired and bio-mass-fued plants use 
supplementaiy firing with gaseous fuel such as methane; others mix coals, 
fi^uently to adjust composite sulllir content; others mix hydrocarbon fiiels such 

30 as coal, shredded tires, trash, and the like. The Input/Loss Method may account 
for all non-coal type hydrocarix)n fuels, but they are assumed to have known 
chemistries with known fuel flows. Mixed fuels (such as multiple coals) are input 
to the system simulator (e.g., EX-FOSS.EXE) as a composite, combining them by 
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vising the FUEL^XE program, although the resultant computed fuel chemistry, 
after proper weighting, is assigned to one fuel (e.g., the most variable coal), by 
assuming the other fuels have fixed or predictable, and known, chemistries. 
Alternatively, such assignment may be made to a combined subset of fuels. 



5 xa, =T^-a^^ (30) 

x(a2 + aj) =T^ (31) 

x(03 + 02/2) =roa (32) 

xa^ =rcoi (33) 

=rso2 (34) 

10 xa.o " - r'Ash - x(a, + 02 + 2aj + a^/l - 1.0) (35) 

where: = 100 - (d^,, + e^c + f + k^„ + l + m/2 + q/2 + 1 + u) 

- m^j,R^ - i.oy[j?^(<|)^ + 1.0)] (36) 

rH2o = (JA«- ^aP) + f+ZPAt/2 + ZP^U/2-bz-bA 

- yi{ZKaJ2 + a,) + bpLjOz (37) 

15 ro2 = d^ + eA«/2 + (G^^-aP) + (l^-Wm 

+ Icac + 3//2 + m^ + p/2 + q + 2r - a - h/l - h^/l 

+ (3o - 2 - 2y + oz)bpLs/2 - x(a7 + ag/2) (38) 

rco2 = dAc. + Cac + YP^t + YPbU + V + w - (1.0 + Y)bpuj 

-xCYRoo + o^ + a,) (39) 

20 rso2=kAe. + / + YbpLs-xa9 (40) 



r'Ash = a-(GAci-ap) + e^2 + fi'2 + w + v-bpLs/2-//2-m/2 
- p/2 - q - 2r + x(ao(YR - ZRy2) + + a/Z + a,/2] 

(41) 

In these relationships the subscript "Act" means an effluent measurement or 
25 assumption(an "actual" value). Theterm JActinEqs.(37)and(38) reiatingtothe 

moles of effluent HjO (or jAct'^Aci referring to air pre-heater inlet conditions). 
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could be iiq>ut as a constant value or measured, given the option suggested by 
error analysis. All other values are either evaluated explicitly based on input data, 
internal models and/or have minor import but carried in the formulations for 
consistency. Eq.(3S) is numerically correct, however if it is reduced by 
S substituting x with and aj with a,^^, see Eqs.(SO) and (5 1 ), it is seen that 
Eq.(3S) offers a trivial solution for aMAF.to; this result illustrates the null effects ash 
has as a neutral substance, and thus its resolution difficulties. 

As a group, these relationships are of critical importance for understanding 
the Input/Loss Method If ftiel chemistry is to be resolved, thus heating value and 

10 thus accurate fuel flow, etc., then stoichiometric relationships generally 
representing Eqs.(30) to (41) must be resolved. These equations are not unique 
in their grouping of terms, further reductions and/or complexities are possible. 
The groupmg of trans adopted here is principally the right-side terms of Eq.(29). 
For the following discussion, assume that fuel ash has zero concentration. 

15 With this assumption, Eqs. (30) through (34) yield five equations with eight 
unknowns. Unknowns include a, through the ienm "a" and ^'x". The term * V 
is a convenience term and could be divided through changing the base of Eq.(29) 
to unity moles of fiiel, thus eliminating use of terms comprising two unknowns, 
x and any ttj. The requirement Ej^..9aj = LOO, assuming ash is temporally 

20 ignored, is again a convenience; again the entire Eq.(29) could be divided through 
by or xa4, setting a base of unity moles of fiiel carbon. Effluent N2 could be 
resolved by difference assuming 100 moles of gaseous effluent (COj, HjO, Oj, 
SO2, and minfM* pollutants are measured or assumed), or could be measured 
directly. By assuming a constant value for fuel nitrogen, a„ a minor fiiel 

25 constituent, theunknownsarereducedto five with five equations. However, close 
examination will quickly indicated that this mathematical system invokes a stiff 
matrix. If 03, agiiin a minor fuel constituent, is assumed constant then Eqs.(3 1) 
and (32) represent two equations with two unknowns, a, and Eq.(33) could 
be solve outright resolving its one unknown These manipulations are 

30 discussed to emphasize that algebraic manipulations must address the physical 
reality of the thermal system as done by the Input/Loss Method. The resolution 
of fuel carbon and fuel hydrogen, of critical importance in resolvmg fuel 
chemistry, should address their intrinsic chemical relationship as a hydrocarbon. 
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and the influoice played by the principal effluent measinmetits^ 4^ J ah • 
Invoking such relationships reduces the Mediod's sensitivity to instrumentation 
errors- The Preferred Embodiment employs this q)proach, relying on die belief 
that coals within a 'generic cat^ory Qlank) and geogrsphical region, have well 

5 defined MAF carbon to MAF hydrogoi relationships. Further, the fuel mole term, 
X, is used in the Prefened Embodiment for this single variable appears in all 
stoichiometric conversions to mass flows (the consistent determination and 
conservation of which is an object of this invention); and represents a iterative 
parameter for Eq.(29). Further, the Preferred Embodiment does not require that 

10 the mmor fuel constituents be assumed constant. Further, as will become 
apparent, the Preferred Embodiment allows use of multidimensional minimization 
techniques which addresses instrumentation errors. 

Returning to the problem as posed by Eq.(29), including fuel ash, if fuel 
chemistry is to be resolved then fundamental problems require solution - 

1 5 independait of algebraic manipulations. These problems include the following: 
1) for eveiy unknown fuel quantity, Oj, there is always, at least, additional 
unknowns in the combustion air term and the fuel moles, "a" and "x"; and 2) the 
fliers solid non-combustible minerals (conunonly referred to as ''ash") may not 
be resolved tiuxiugh effluent measurements since it is a neutral non-gaseous 

20 substance, the same inpvt of xaio in the fuel appears in the effluent, further, its 
outiet flow from the system is not possible to measure directiy witii any reascmable 
accura^. 

The first of these problems is solved by reducing quantities to a molar 
MAF bases, eliminating flie influence of the two components not chemically 

25 involved in the hydrocarbon fuel - watCT and ash • and then solving for Omaf^ 
terms algebraically employing correlations of molar MAF fractions of fuel 
hydrogen, nitrogen and oxygen as a function of carbon. Fuel MAF carbon is 
resolved explicitiy by first solving Eq.(33) for x, converting to a MAF base, 
subtracting Eq.(32) from Eq.(31) to eliminate the XMAFaMAF.2 t^s, and into the 

30 result substituting Eq.(33ys solution for Xmaf* yielding: 

aMAi.4 =-rco2(A5-2A3)/[rco2(B5-2B3)-r«o + 2roJ (42) 
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Note that the w & V terms of Eq.(39X related to reject & refuse carbon may be 
expanded, as functions of ocmam (see Eq.(29) notes), if their magnitudes afTect the 
accuracy of Omaf^^ Given Omam» ^ resolved via Eq. (56), based on Eq. (33). 
In Eq.(42) the constants and Bj are developed from reference fuel data 
5 associated with coal samples taken fiom the actual As-Fired, or obtained from 
generic specification of the coal Rank: 

Omaf-i ==A,+B|aMAM ^^^y 
Omaf-j ^^Aj + BjOmaf^ (^) 

10 Fuel MAF sulfur is resolved e3q)licitly by solving Eq.(34): 

<^hAA?'6 "1^802 '>^IAF ^^^^ 

Fuel MAF water is resolved explicitly by solving Eq.(31) for x, and substituting 
into Eq.(32), yielding: 



ttMAF.2 (^MAF-jT^mO • 

aMAF-5ro2)/(ro:-r««/2) (47) 

IS Both fuel MAF nitrogen and OT^gen concentrations are small, and 

typically may be fixed as constants, or described by Eqs.(43) and (44). The MAF 
sulfur concentratim may also be computed assuming aMAF-6 ^ A^ ^ B^ocmaf^, if the 
effluent SO^ measuremmt is questionable (or if Umestone is injected and its 
conversion rate is questionable). As a further altomative, either the o>g^gen or 

20 sulfur could be used to assure that Tr^^^j^ jx^^^ = 1 .0 (assuming a^AP^ "maf-?* 
^maf^ «mafj9i any other secondary fiiel has zero or known concentration). The 
sole criteria in deciding the exact methodology is the reliability and availability 
of efiOuent data and its relative impact on MAF fuel terms. For example, nitrogen 
is a major efQuent, however using Eq.(30) to solve for fuel nitrogen, typically a 

25 vety minor component of the fiiel's mak«ip, would invite even slight errors made 
in determining eflQuait Nj (made either by direct measurement or determined by 
difference suggested in Eq.(36)) to greatly amplify the imcertainty in Omaf-i- 
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Indeed, the Fiq term is used to resolve the combustion O,, the term ""a** via Eq. 
(58); wherein the term Xmaf«maw is resolved toough Eq^43)such tbdX Tfo, wWdi 
is dependent on ''a**, may be resolved Besides carbon, an important tenn is the 
fuel MAF hydrogm concmtration which must represent a valid functionality with 
5 carbon. Experience has shown valid functionality may be easily adiieved since 
the hydrogen:caibon chemical bond is predominate with carbonrcarbon; indeed 
the variability of a specific coal Rank lies not in its base chemistry but rather in 
its water and ash. FIG. 3 is a sample plot of MAF hydrogen versus MAF carbon 
for a coal fuel having highly variable As-Fired heating value, water and ash 

10 contents. In FIG. 3, both axes are MAF molar fractions per mole of MAF fueL 
To this point in the process, resolution of fuel chemistry has been possible 
by assuming a Moisture-Ash-Free (MAF) basis. In carrying this out, it is not 
critical to the invention that carbon be computed as the independent quantity. 
HydrogOT, or any oAer MAF fuel element could be considered as the independent 

15 variable. For example, if solving for hydrogen first, then Omata =/(aMAF.5)- 
However, the preferred procedure is to place indq)endence witti carbon, v^ere of 
course the greatest accuracy and sensitivity may be fowd in the efQuent 
measurements. CO2 does not exist in the combustion air to any appreciable 
concentration, it does not leak into the system, it is generated only from 

20 combustion and as a major effluoit has obvious sensitivity to resolution of aMAF-4 
ofEq.(42). 

The second problm, dealing vdth ash, requires examination of the total 
system. The only system effect of fuel ash is as a pure dilutive or concentrative 
influence on fiiel, and of course on the heating value. From a qualitative 

25 viewpoint, as fuel ash increases at the expense of carbon (for example), the 
amount of combustion air required to produce the same effluent Oj, given less 
carbon to combust, will deo-ease. Given a decreasing heating value (higher ash) 
increased fuel flow is required to meet the same useful energy flow. Thusanideal 
system parameter for such sensitivities, which is routinely measured at fossil- 

30 fueled systems, is the Air/Fuel ratio. Generally such sensitivities are reasonable, 
a 10 percent increase in ash for a common coal will cause a 5 percent decrease in 
the Air/Fuel ratio, although many system and fuel parameters may affect such 
sensitivities. The wet, mass base, combustion Air/Fuel ratio (termed AF), as 
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computed by the system simulator is developed as follows: 

^=(niAir+nWBm)/ni^p (48A) 
^F=(l+P)[(a + a<|)^NA^+bAN„«J/(xNAF) (48B) 

Expanding the tenn xNaf in Eq.(48B): 

5 xNAF=x(Z;^^jOj + N,oa,o) (49) 

and then employing the following definitions of MAF fuel moles and fuel 
constituents: 

Xmaf = X / (1 .0 + aMAF-2 + ttMAF-Io) (50) 

Omaf-j = aj(1.0 + cwF.j+aMAF-ia) (51) 

10 allows substitution of Eqs.(SO) and (5 1 ) into Eq.(49) for x and Oj, cancelling the 
(1-0 t^^MAFo ttMAF-io) as intended, and then substituting into Eq. (48B) 
yields a solvable form: 

XNaf = XMAl<J^jH)...9NjaMAF^ + NioOmawo) (52) 
AF = (1.0 + P)[(a + «|>>JNa, + KVlnol 
15 /[W5W.9Nj«MAi4+N,oa„AF.ia)] (53) 

Simplifying Eq.(53) and solving for MAF fuel ash. Omaf-ioi yields the following 
results. Note in Eq.(54) that a normalized Air/Fuel ratio is used, noraialized to 
actual plant data, defined by Eq.(57). Xmaf is substituted using Eq.(56). 

20 where: s (i.o + P)[(a + a<l>Aa)NAi, + bANH2Ql/^A« (55) 

^MAF ~ ^C02 ^ *MAF.4 (56) 

AF^ s AF^ (AF^/AF^ (57) 

a =(rNj-XMAFaMAF-l)/<l>A« (58) 
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The variable AF^^ is the wet Air/Fuel ratio from the system's data collection 
devices (an unsealed signal); AF^ is the wet Air/Fuel ratio consistent with 
reference fuel chemistry; and AF^^ is the w^ Air/Fuel ratio consistent with 
reference fuel chemistiy but as would be produced torn the system's data 
5 collection devices (a signal used for scaling). The value of N,o in Eq. (54) is input 
as a constant, or fitted as a function of Omaf-io (*us solving a quadratic equation), 
or fitted as a function of HHVm^. Note that a system's indicated fuel flow 
measurement could obWously be used in place of ^F, 2q>plying similar techniques 
as demonstrated in detennining Omaf-io. However, use ofy4F is preferred since it 

10 integrally involves effluent and combustion air terras (through and T^, 

and thus through such dependencies allows error analysis techniques to be 
operational and practical. 

It is noteworthy that the explicit procedure of determining fuel ash, and 
tfirough use of die term (1.0 + aMAF.2 + aMAwo) of Eqs,(50) and (51), allows any 

15 errors made in fiiel water, Omaf-ji to be off-set by fuel ash, Omaf-io This must 
occur since any given quantity xa| (wet-base) must be equivalent to Xf^^uAv^ 
(MAF-base); if not, such wet to \fAF conversions would numerically cause 
inconsistmcies m the conqsuted Air/Fuel ratio. There may be only one AF - wet- 
base is the same as MAF. Furtho*, the effect may easily be demoiistrated through 

20 off-line study. 

The selection of specific Calculational options is depmdent upon the 
characteristics of fossil fuels. As an example using coals of various Ranks and 
other fuels, FIG. 4 illustrates a table of options relating to fiiel water, fuel ash and 
the fuel's typical variability. Note that these options are subject to application of 

25 the multidimensional minimization techniques described herein; for example an 
assumption of constant effluent HjO may be corrected using error analysis. 

Thus all fuel constituents are therefore determined on a MAF bases. Given 
these values, the wet base fiiel molar fi^ons are determined, the wet base moles 
of fuel (x), followed by wet base (As-Fired) weight firactions: 



ttj =aMAF^/0.0 + aMAF.2 + aMAF.io) (59A) 

X =XMAF(iO + aMAF.2+aMAF.lQ) (59B) 

WF^ ^aiti^/CDxfi) (60A) 

»^<iry-i = »W^i/(1.0-»7y (60B) 
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The present Input/Loss M^od has important advantages over the prior 

techniques discussed in tfie Badcground. The foiioi^ong table clarifies differences. 

Determinations of 
Fossil Fuel Chemistries 



Method 


Water 


Ash 


MAF Elements Comments 


Input*Output 


no 


no 


no 


Chemistry not mentioned. 


Heat-Loss 


no 


no 


no 


Chemistry assimied 


Output/Loss 


no 


no 


no 


"Approximate" chemistry assumed. 


Buna 


no 


no 


no 


Chmistry assumed. 


Fuel Water & 










Ash Meters 


yes 


yes 


no 


MAF chemistiy assumed known, 










may adjust MAF heating value. 


Thermal 








Analysis 


yes 


yes 


yes 


Off-Line, laboratory use with high 
accuracy gas measurements^ 
with »gram samples. 


Munukuda 


yes 


no 


yes 


No air leakage, fuel ash is assumed 
known. 


Pat 5367470 


yes 


no 


no 


Water iterated, known diy 
chemistry relationships. 


Pat 5790420 


yes 


no 


no 


Water iterated, known dry 


Input/Loss 








chemistry relationships. 


yes 


yes 


yes 


Explicit calculation for water, ash 



and carbon, vntti error analysis. 

25 Heating Valyg Calculations 

Variations to a refenaice heating value may be determined by calculating 
the diflferential effects of MAF reference versus actual MAF fuel chemistry. The 
As-Fired fuel chemistiy could be used to produce the MAF or dry. The heating 
value calculation involves use of a correlation relating oxygen, carbon, hydrogen 

3 0 and sulfur contoits to a MAF heating value, correcting for ash to obtain a dry-base 
heating value, then finally correcting for water to obtain the As-Fired value. The 
As-Fired value is then corrected for a constant pressure process if qjpropriate. 
The preferred correlation is sciwitifically based on the chemical binding energies 
between hydrocarbon elements. The correlation is taken from the worics of 

35 Ghamarian & Cambel, which is based in-part on the well known work of Szargut 
and Szargut & Stryrylska. The references include: A. Ghamarian & A. B, 
Cambel, "Energy/Exergy Analysis of Pressurized Fluidized Bed Combustor'', 
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Proceedings of the Intersociety ^eigy Conversion Engineering Conference, 
August 8-12, 1982, pp. 323-327; A. Gbamarian & A B. Cambel, "Exeigy 
Analysis of Illinois No. 6 Coal", Energy, Vol. 7, No. 6, 1982, pp. 483-488; J. 
Szargut, "Intematiottal Progress in Second Law Analysis", Energy^ Vol. 5, 1980, 
5 pp. 709-718; and J. Szargut & T. Stryiylska, "Approximate Determination of the 
Exergy of Fuels", Broinstoff-Warme-krafi, Vol. 16, No. 12, December 1964, pp. 
589-596. The correlation is accurate to within ±0.7 percent AHHV^iAF-Ref 
deviation for over four dozen short- and long-chained hydrocarbon compounds. 
For hydrocarbons like coal, demonstrated below, normally having a low oxygen 

10 content, the correlation's accuracy is estimated at ±0.5 percent A similar 
correlation by the same authors exists for hydrocarbons wth high oxygen content 
To further reduce these errors, the Input/Loss Method computes a term dHHV^. 
detai based on a referraice MAF-base heating value, HHVMAjLRrf and the reference 
concentrations of oxygen, carbon, hydrogen and sulfur (established during 

15 initialization). The term HHVm^.^ is computed via Eq.(62) based on actual 
MAF cmcadrations of Eqs.(42) thru (46), or otherwise determined as discussed. 
Using Eqs.(63), (64) and (65) the on-line As-Fired heating value, HHVaf, is 
resultant The term Nmaf is the molecular weight of ttie MAF-base fuel (without 
watfflr and ash). 

20 AHHVMAM.h.=HHVM«,^^ 

-(- 178387.18aMAF.3+ 183591.92o„^^ + 78143.68aMAF-s 

+ 127692.00aMAF^Rrf/N»^^^ (61) 
HHVma^»c«i. = (- 1 78387. 18aM^.3+ 183591.92oMAM + 78143.68a,^F.5 

+ 127692.00Omaf.6)ac.u,i/Nmai^a«o., (62) 

25 HHVmaf = HHVM,u:^+AHHVMAMdu (63) 

HHV,^ = mW^l.O -WF^t^ (64) 

HHVaf = KHVa^LO-fFT-j) (65) 

The preferred correlations used to determine heating values for the present 
invention are based on chemical binding eneigies. Studies have demonstrated that 
30 traditional correlations, such as the Mott-Spooner correlation based on Dulong's 
formula - well known in the industry - are not adeqiiate. The Preferred 
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Embodiment of the present invention requires at least die coefficients used in 
determining heating value to fall within cotain ranges associated with three 
principal constituents of coal. Studies have indicated that using the above 
prefeired constants, which fall within die required ranges, reduces the standard 

5 deviation of five dozen wildly varying coal analyses bom ±228 todb92 ABtu/lbm. 
The ranges of these coefficients, i.e., multiples the molar firactions ttj in Eqs.(61) 
and (62), or their equivalent wei^t ftactions (for this presentation of ranges, die 
symbol WF^ represents percent weight of j), include the following: for carbon 
molar fraction I74160ac«rton/Nfl.ei lM970a„rt«/Nfiieb or in weight percent 

10 carbon, MSJFFe^to 154^^;^; for hydrogen molar fraction 776 lOahy^fa^^j^ 
to 114910ahyd„,gJNfi,ei assuming molecular hydrogen, or in weight percent 
hydrogen, 3%5fVF,,^^ to SIOWF,,^, and for oxygen molar fraction- 
163 190ao,yg^fi^ to ' 1 82390a,,,yg^f^ assuming molecular oxygen, or in weight 
percent oxygen, -5 1 WF,^ to -57 WF^^^. These ranges are independent of the 

15 fiiel base, whether MAF, dry or As-Fired fuel constituents are used. 

Thus the As-Fired fiiel heating value has been determined. If a correction 
for a constant pressure process (HHVP) is required to the computed HHV^p, this 
correction may be made U3ing the procedure of PTC 4.1. Iterations of the above 
system, combustion, thermal efficiency, fuel chemistry and heating value 

20 calculations are preformed repeatability converging fuel moles, x, the higher 
heating value, HHV^, such that total fuel energy flow, As-Fired, is adequate to 
meet the Useful Energy Flow Delivered, and minor stoichiometric terms to within 
tight tolerances. 

Gross Error Detection and Mas s Balances 

25 After convergence, the Input/Lx>ss Method employs the primary 

measurements d^^,, and AF^ (used to determine fiiel dimistiy, thus heating 
value, etc.) in comparisons to the computed values, Jqj, and AFq^ determined 
&om the system simulator, to deduce inconsistencies in measurements. However, 
it has been found that such comparisons indicate only gross inconsistencies, and 

30 not necessarily those associated with incorrect system parameters such as air pre- 
heater leakage (R^ or combustion air Oj (A^J, and when used in an on-lme 
application were not of value. A more robust and operational error analysis 
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employs multidimensioiial minimization techniques, discussed below. 

It is an important objective of ttiis invention to assure conservation of mass 
flows. Such conservation, using the terms developed by this invention, 
int«x:onnecting boiler efficimcy, shaft power and useful energy flow and 
S effluents, assures the thermal system operator that consistency of these 
computations is acceptable. Mass balances are specified in the following table, 
note that if inlet and outlet mass flows disagree by more than 0.2%, errors are 
considered significant. Such errors may be associated with computer 
programming mistakes and, most importantly, inconsistencies in the principle . 
10 effluent measurements. 



Fuel Flow Rate, m^p, Eq.(21) = BBTC/[nboii«(HHVP + HBC)] 
Combustion Dry Air Flow Rate = mAF(l+P)(a + ac^Ac^N^iAxNAp) 
Combustion Air Moisture Flow Rate= m^j^l + P)bAN„2o^(xNAF) 
In-Leakage of Water and Steam = m^fbJtiwJO^AF) 
15 Pure Limestone (CaCOj) Flow =-ja^\ + Y)lW^3MAia-- 

EINLET FLOWS 

Dry Gas Flow finom the Combustion 
Region (i.e.. Boiler Effluent, 
assuming 100 mole base) - mf^lOOHoJiRj^^ xN^) 
20 Dry Air Leakage Flow at Boundaiy = m^pPCa + ^^Aj^Ai/O^AF) 

Combustion Moisture plus Air Leakage 

Moisture at Boundary = TOafO PbA)NH2o^(xNAF) 
Reject & Refuse Carbon Flow = mAi<v + w)Nc/(xNaf) 
Calcium Sulfate with Wat^ Flow = inAi^bpLsNciso4.2H20^(xAFNAF) 
25 Calcium Oxide Flow = m^l - a + y)h^aic^iiiNj^) 

Ash (Bottom, Fly & Dust) Flow ^jnAFBioMUsh^F 

EOUTLET FLOWS 

EfflwgntCalwlations 

By knowing the fuel chemistry and flow rate, and the complete 
30 stoichiometric relationships, calculating individual emission flow rates, mspecies-i 
(ib/hr), may occur as follows: 
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in,pec«.i=niAF*iNi/[xNAF] (69) 

where *i is the molar firaction of an eflBuent species on a dry-basis. The tenn *j 
derives directly from solutions or measurements of the right-hand terms of 
Eq.(29), for example Oso2 Thetermxisthemolesof fiiel permoleof dry 
5 gaseous eflQuent The emission rate per species, in units of pounds per million Btu 
of fuel energy input, termed ERj, is given by the following: 

ERi = lO^m^ieJCmApHHY^F) (70A) 
= 10^OiN/[xN,,pHHV^] (70B) 

The emissions rate may be evaluated independently of the As-Fired total fiiei flow 
10 rate, Eq.(70B). However, the computational accuracy of the fuel flow rate, as 
determined using the present approach, intrinsically affects the emissions rate 
through HHV^ x and N^p- Further, the process described herein allows the 
determination of total dry volumetric flow, at standard conditions of gaseous 
effluent, denoted by VF, as required by environmental r^;ulations. VF is 
1 5 determined by the following (in standard-fl^/hr): 

VF^P^A^JlxHj,,] (71) 

where Pga, and Ng^, are the standard dwisity and average molecular weight of the 
effluent dry gas. Of course, to d^nnine the mass flow of all effluents Eq.(69) 
may be summed 

20 Minimization of Instrume ntation RrmrR 

If the aforementioned mor diecks exceed an error tolerance, or for routine 
analysis of instrumentation errors, the Method of this invention may anploy 
multidimensional minimization tedmiques to adjust measured and/or assumed 
effluent data (d^ & J^, and measiired and/or assumed system data (i.e., AF^ 

25 Racv <l>Act). resulting in revised fiiel constituents. Although in the system of 
equations of the Prefmed Bnbodiment, the effluent is considered a known and 
an accurate measuranent, having rK> discemable error, this is not requii^ as wUl 
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be observed in the following discussion. The assumption of known efQuent O2 
is supported inpractice since industrial O2 instrumentation, commonly employing 
zirconium oxide technology, is relatively inexpensive. Thus a given thermal 
system may install numerous O2 instmments, upstream and downstream of the air 
S pre-heater, to achieve a result free of discemable error. Furthermore, the & 
<|)Act parameters, important components of the combustion equation, may also be 
determined using O2 sensors. 

When using the techniques of the Preferred Embodiment for minimizing 
error in the air pre-heater leakage factor Oj measurements are required 

10 between the heat exchangers/combustion region and the air pre-heater, see Fig.6, 
and preferably immediately downstream from any effluent ash removal equipment 
but before the air pre-heater. 

After converging on fuel moles and heating value, with a given effluent Oj, 
the Method of this invention then utilizes certain parameters which are unique to 

15 a given fuel to determine instrumentation errors. These parameters include the 
moisture-ash-free (MAP) heating value, and certain parameters (termed L^) 
developed from Roughton's work of 1980. Study of large coal data bases has 
indicated a remarkable consistency in MAF heating value for a given coal of the 
same Rank and geographical region; for example, a Pennsylvania sub-bituminous 

20 coal's HHVmaf P^r cent deviation was found to be less than 0.02%. Roughton 
found that the 1^^^ parameter was within ±1 .0% with 95% conOdence for "a wide 
range of powar station coals". This has been generally confirmed, but again 
confined to individual coal Ranks and regions; variations of ±0.5% for unique 
coals have been found to be common. The parameters for the total fuel, fuel 

25 water and fuel ash are determined in the following manner: 

• ^0m^H2O • ^WF-theor^MAF-loNAshl ^ (Xory^heorNory-FuclHHVD^) 

(72) 

Note that: ^^0,y.^Diy-fvJ^Oiy= >WF^heofNMAF-FMriHHVM;^. 



30 W = Jcheo^HlO / (XD^.he«ND^ue^fflVo^) (73) 
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where J^^^ k^, and are used to develop these terms as detennined 
directfy fixim die iq>dated fiiel chonisto^, assun^ Note 
that no air pre-heater leakage is assumed for theoretical combustion: J? = 1.0, 
p = 0.0, thus J^^ = j^^ and = 0.20948. 
5 The prefinred procedure then establishes correlations of the Zr, parameters 

as functions of the actual input effluent or system data. For example, the 
following correlation relates the computed fuel factor Ijraei to effluent CO^: 

= K„ + Kiid^ct + K,3d^ (75) 

Eq.(75) is established by varying COj (dj^c) and calculating Lf^ for three cases, 
10 keeping all other effluent parameters constant. Such variation is done 
automatically so that the reference value, L^.^^ is encompassed within the 
correlation's range. K,„ K,2 & K,, are associated fitting constants. In like 
manner, correlations are formed for other parameters: 
(Lw«er+^Ash) versus J^tt 

15 (Avaier + ^Aah) VCrSUSiiF^^ 

versus air pre-heater leakage factor (the term /?Act) 
Lpuei versus concentration in the 

combustion air (A^, vrtiere A^ =1/(1+ 

The Ii4)ut/Lx)ss Method determines all such correlations automatically including 
20 their proper ranges. Note that the parameter L'^^^ was found to yield best 
sensitivities to effluent COj, as opposed to Roughtorfs evaluation of the same 
term, \Axich is used for sensitivity with and A^^. He defined Z^,, based 
on dried fuel, as the mass flow of diy gas produced firom theoretical combustion, 
Xi>yW^D«y^ divided by the total energy of the dry fuel, Xo^^hecNo^^fflVD^. 
25 For perfect effluent measurrnirats, without error, this later form and that used in 
Eq.(72) are of course identical, they represent different sides of the same 
theor^cal combustion equation using dry fuel. Use ofl^^^ and L'^^ was found 
of importance to uniquely sensitize their respective terms. The Preferred 
Embodiment of this invention employs the summed parameters (L^^ + as 
30 afimctionofJActand^F^individuaUy;vereusIw«er=/(JAJ A^d> =/(^/Vvct)- 
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This summation was found essmtially constant for a given fuel, more so than its 
individual terms. The Preferred Embodiment of this invention also miploys, in 
addition to the Zrj parameters, the filers MAF heating value. Again, the MAF 
heating value has obvious hnport given its constant nature for a given coal, but use 
5 of the diy and/or the As-Fired (wet) heating value could be applicable if sample 
coal was especially prepared for testing purposes leading to correction of 
instrumentation signals (corrections being applicable to normal operations using 
fuel with unknown chemistry and heating value). 

These terms then form individual error functions, Fj, the minimization of 
10 their aggregate forms the Objective Function, solved using well known 
techniques. 

^i = [(A-iRcf.iy^iufJ'-[l-0-nn-iAi(^ (76) 
Objective Function = E Fj (77) 

where are computed values based on actual data, are corresponding 

IS reference values; and HHV„ are computed heating values, HHVg^^ are 
corresponding refermce values for MAF, diy and/or the As-Fired (i.e., n = 1 ,2 or 
3). With tfiis formulation, the power plant operator has the option of choosing any 
one or all of the five efiluent and system measurmients, with up to three types of 
heating values each. Further, any other relevant effluent and/or system 

20 measurement could also be incorporated, such as effluent Oj [L^^ =/(gAct)]» 
effluent SO2 [Z^ =^f(KJ)l and/or indicated fuel mass flow [L^ ^/(m^?^^]. 
It should be noted that this error analysis would not be possible without the EX- 
FOSS program's ability to compute consistmt effluent data based on given inputs 
of fuel chraiistry and system parameters. Measured efflumts are not input to EX- 

25 FOSS, but orJy effect computed fiiel chraiistry. Capitalizing on this, the Method 
of tiiis invention then forces consistency with fuel chemistry by correcting 
effluents (and thus chemistry) and system param^ers through error analysis. The 
end result is consistent fiiel flow, effluent flows, emission rates, boiler efflciency 
and system efficiency. 

30 Note that the term ^ference fuel characteristics** includes an average or 

typical fuel chemistry and associated MAF heating value, preferably based on 
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historical data collection of ultimate analyses. The term ''elementary" when used 
herein, referes to a fuel's composition by individual elemmts (denoted by the term 

ajinEq.(29)). The values of £,^^^,XKef.watBr^d£Kef.Asii^^ ^^'^ 
of the reference fuel charactoistics, computed using the reference fuel chemistry. 
S Reference fuel characteristics also includes Aether the variability of fu^ 

and fiiel ash in the As-Fired condition is predictable, or not For any given fuel: 
fiiel water may be held constant (including zero); ftiel ash may be held constant 
(including zero); an accurate functionality may be observed for either or both 
[e.g., a^i^io = /(HHVm^)]; and/or fuel water and/or fiiel ash may be treated as 

10 unknowns, determined using procedures of this invention. All of these possible 
variations for the treatment of fuel water and ash are included as a portion of the 
Method. Reference fuel characteristics also contain reasonability limits of the 
computed elementary fuel constituents. Reference fiiel diaracteristics also contain 
fitting constants (e.g., Ki,, & K^) associated with ail correlations relating 

1 S dependent quantities to system measurements. 

It was found that the Preferred Embodiment used for the minimization of 
the Objective Function is the Conjugate Gradient metlKKl, described in the book: 
W.H. Press, S.A. Teukolsky, W.T. Vettering & B.P. Flannery, Numerical Recipes 
in FORTRAN 77. The Art of Scientific Computing. Cambridge University Press, 

20 Cambridge and New York ( 1 992). This technique is sufficiently robust as not to 
require side constraints such as limiting boundary conditions. It produces effluent 
and system parameters (e.g., COj^ HjO, R, etc.) which as a set of data represents 
mmimum cnor, thus correction factors are detmnined: for the CO2 signal 
dcc/dAo, for the indicated Air/Fuel ratio AF^^ = AF^^ fAF^ etc. This 

25 technique, as others, require die partial derivatives of Eq.(76), thus: dFJdd^ 
3FJdJ^ dFf/dAF^ dFJdR^ and dF/dA^ are computed in a routine manner. 
Other published techniques are common, such as those described in the book: 
G.N. Vanderplaats, NumCTical Optimization Te chniques for Engineering Design. 
McGraw.HillBookCompany,NewYoric(1984). Given the MetfKxi's procedure 

30 of allowiqg for effluent and system measurements to be compared to 
corresponding, but consistently computed values based on updated fiiel 
constituents, common minimization techniques involving constrained searches, 
could be alternatively applied. These include Sequential Linear Programming, 
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but based on the computed values dcai, Jen and ^ca- Obviously the deviations 
between the measured and calculated values of Omaf^ therefore m^ be 
determined: 

= (ttMAF^E - ttMAF/aw) I «MAF/ai* (78) 

5 Any non-zero value of Eq.(78) may be ultimately attributed to inconsistencies in 
the effluMit and/or system measurements. Therefore a valid fuel composition is 
the one which minimizes these deviations thus being consistent with both 
stoichiomeby and measur^ents. This is not preferred since no use is made of the 
/?Act ^d terms as might affect effluent measurements and/or independent 
10 parameters (e.g., L^^ & HHVm^). 

Summary of Thermodvnamic Understanding 

The present Input/Loss Method provides imjroved understanding of the 
thermodynamics of fossil systems, as compared with the prior approaches 
discussed in the Background. The present Input/Loss approach allows 
IS determination of complete As-Fired fuel chemistries, heating values, fuel flow, 
effluent flow, consistently calculated ash flow, consistmtly calculated emission 
rates, and consistently calculated system efficimcy. Comparisons to other m^ods 
are presented in the following table: 

Complete As-Fired Calc. Calc. Calc. 

20 As-Fired Heating Fuel Effluent Ash Emission System 

Mfilllfid FuglChem, ValW Flow Flow Flnw Ratef; Rfficiencv 

Input-Ou^ut NA PM PM NP NP NP yes 

Heat-Loss PM PM NA NP MP NA yes 

Output/Loss (no) PM yes PM NP yes yes 

25 Buna PM PM yes NA NA NA yes 
Fuel Water & 

Ash Meters (no) (no) NA NA NA NA NA 
Thermal 

Analysis (no) (no) NA NA NA NA NA 

30 Munukutla (no) (no) (no) NA NA NA yes 

Pat. 5367470 (no) (no) (no) (no) NA yes yes 

Pat 5790420 (no) (no) (no) (no) NA yes yes 

Input/Loss yes yes yes yes yes yes yes 
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Notes include: 

PM => Perfect measurement or perfect assiunptions made. 
NA => Not Applicable, or not done given original presentation 
of method. 

5 NP => Not Possible. 

"Calc " => Calculated with consistency based on non-^ect 
measurement or assumptions, 
(no) -> Results are possible but restricted to either known MAF fuel 
chemistry, and/or known fuel ash, and/or known fuel 
10 water; for Munukutla, his method is restricted to 

constant and perfect assumptions as to fiiel ash. 



At the practical level of applying the Input/Loss Method at an operating 
power plant» its use of an air pre-heater leakage factor (/?Act)» ^n air pre-heater 
dilution &ctor (p), and a variable Oj/cmnbustion air term {Aj^ thus 4>^, which 

15 allow use of effluent measurements on either side of the air pre-heater, coupled 
with techniques which correct efQuent instrumentation through error analysis, 
assure operational and reliable thennal performance monitoring. Considering that 
boiler efficiency is formed from fiiel constituents (based on measured effluents), 
that a set of consistent effluents are computed based on these fiiel constituents, 

20 that the measured effluents are corrected to yield minimum differences with the 
computed through error analysis, the power plant operator has assurance of 
consistency. Thus this Method, considering its procedures, defines "^self 
consistency**. 

When q>plying the MetfKx] of this invention for the inqirovement of 
25 thermal system efficiency the Preferred Embodiment is through an established 
technique using Fuel Consumption Indices (FCI). FCl are developed from 
exergy and theimodynamic irreversibility computations of the system's major 
ccmiponents and processes based on the Second Law of thermocfynamics. These 
components and processes include in part: the combustion process, combustion 
30 gas to woiking fluid heat exchangers, the turbine cycle, and the process of direct 
generation of electricity . FCls are computed by the EX-FOSS computer program. 
For the purposes of this invention this technique is defined as a ""Second Law 
analysis of the thermal system". The most important inputs to this technique 
includes specification of As-Fired fuel chemistry, As-Fired fuel mass flow and 
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combustion gas mass flows; which given the consistent approach of the present 
invention, are fully resolved Other required inputs include system mass flows 
including those of the working fluid, system enthalpy flow balances and 
miscellaneous data. Thus the Method of this invention are uniquely suited to 
5 allow the benefits of Secrad Law analysis of thermal systems to be implemented 
The reference for this technique is: F.D. Lang and KLF. Hom, 'Tuel Consumption 
Index for Proper Monitoring of Power Plants", EPRI 1991 Heat Rate 
Improvement Conference, May 7-9, Scottsdale, AZ, also available fi-om Exergetic 
Systems, Inc- of San Rafael, CA in its latest version. 

10 Summary of Procedure 

In summary, details of the procedure involves first several initializing 
procedures, followed by sequential procedures preformed by a computer operatmg 
in communication with the systems' data collection devices v^ile it is operating, 
in real time. 

IS Prior to real time operation, the Iiq)ut/Loss Method requires certain 

initializing data involving refermce fuel chemistry and heating value, and 
stoichiometric data associated with the reference fiiel chemistry. To supply such 
data a series of reference fuel chenustries, as ultimate analyses (identifying by 
weight fractions the mdividual elements), established from historical laboratory 

20 reports is desirable to determine relationships between an independently 
determined MAF molar carbon fractim ct^^ and dependent hydrogen, nitrogen, 
oxygen and sulfur fractions; i.e., aMAF^==/(aMAF^)- In coal fuel, nitrogen, oxygen 
and sulfrir compon^its are small relative to carbon, hydrogen, water and ash, and 
thus have minor unportance to the fuel's heating value. However, the relationship 

25 =/(aMAF-4) should be established based on actual data. An associated 

reference hig^^ heating value, consistently established from the same historical 
laboratory reports acquired for the fuel, is desirable, reduced to a MAF base. 
However, if no such data are available, then the fuel's minor components of 
nitrogm, oxygen and/or sulfur may be held constant and/or calculated explicitly 

30 if effluent data has the appropriate accuracy, while the relationship ocmaf.5 

(^MAF-4)> ^d the reference MAF heating value may be established from coal 
analyses associated with the generic coal type. During the initialization, an error 
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analysis is preferable using sampled data from the system to determine both the 
interaction of instrumentation errors on effluent CO2, effluent HjO, Air/F uel ratio, 
the term (and thus air pre-heater leakage), and the term A, and to confirm 
model selections related to fuel chemistry computations. 
5 FIG« 1 illustrates the initialization processes. Box 1 10 indicates the initial 

(and one-time) setups for the indicated items. Box 120 depicts the process of 
developing a reference fuel, involving all common fuels used by the system. 
Reference data is input as wet-base weight fractions to a system simulator such as 
EX-FOSS JBXE, but also reduced to Moisture-Ash-Free (MAF) base for use by 

10 Eq.(6 1 ). Box 1 30 refnesents a simple organization of all initialization data. Box 
140 depicts using the same system simulator used when on-line, for off-line 
sensitivity studies resulting in selection of appropriate calcuiational options. 
FIG. 3 presents a typical relationship between MAF molar hydrogen and MAF 
molar carbon. This particular data is based on wildly varying coal fiiel, whose ash 

15 varied by ±19 percent about the mean, whose fuel water varied by ±40 percent, 
and whose heating value typically varied by ±12 percent All seccmdaiy fuels 
(e.g., coals and/or gas and/or oil fiiels) are specified, and assumed to have constant 
values. FIG. 4 suggests Calcuiational options associated with various fuels. Each 
of these options should be confirmed durii^ the initialization process involving 

20 error analysis. Box 140 of FIG. 1 . With these initialization steps completed, the 
analysis proceeds to real time, on-line analysis. 

Asusedherein, the tenn "obtained" or ''obtaining** is defined as measuring, 
calculating, assuming a non-zero value, assuming a zero value, estimating, 
gathering from a published reference work, gathering fix)m a data base, or any 

25 other operational approach 

This paragrq>h discusses and defines ""opmting parameters'" obtained 
fit>m a thermal system applicaUe to tfus inventioa . EfQuent CO2, Oj, NOx> and 
SO2 measurements are made (at the systrai's effluent stream or "smoke stack**), 
or bef<m the mr heater (termed the "boiler" side of the air pre-heater, i.e., 

30 upstream of the air pre-heater device). O2 measurements are preferred before the 
air pre-heater. Indicated wet-base combustion Air/Fuel mass ratio and effluent 
H2O measuremoits are required, or assumptions made, dq^ending on the reference 
fuel characteristics and appropriate error analyses. Measurements con4)rising the 
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Air/Fuel ratio could be made on a volume base» or a dry-base, then converted. 
Effluent temperature is measiued, that is the meiBgQ temperature associated with 
tfie combustion gases at the boundary of the system (caution must be exercised in 
measuring non-stratified gas flows). The inlet/oudet ratio of CO2 (preferred), CO, 
S or O2 aooss the air pre-heato* (these could be obtained off-line, based on periodic 
testing or judgement, but on-line is preferredX is used for die determination of air 
pre-heater leakage in accordance with ASME PTC 4.3, using the term (and 
thus P, see notes associated with Eq. (29)). Fuel temperature at an appropriate 
thermodynamic boundary of the ^stem is obtained Air psychrometrics are 

1 0 obtained at the boundary of tfie system (e,g., diy and wet bulb temperatures, or dry 
bulb and relative humidity). The discharge temperatures of tfie air as it exits each 
air heating or cooling device (but before it reacts with the fuel) are required; for 
example, such devices might include the air pre-heater, forced-draft fan, steant-to- 
air heater, etc. Measurements are required to determine the total enei:gy flow 

IS deposition to the working fluid from the combustion gases. For a power plant, 
such measurements typically include feedwater flow to the steam generator, 
feedwater pressure and temperature, steam flow from flie steam generator if 
different than the feedwater flow, steam pressure, steam temperature or quality (or 
assumed quality), and, if ^plicable, reheat flows, and reheat inlet and outlet 

20 pressures and temperatures. For a conventional power plant, determination of 
accurate reheat flows generally requires understanding of steam turbine flow 
distributions (involving shaft seals, feedwater heaters, bypass leakages, 
attemperation spray flows and the like). 

FIG. 2 presents the Calculational sequences associated with the Method 

25 when monitoring a fossil-fired system on-line, i.e., in essentially real time. Box 
2S0 represents data passed from the initialization process dq>icted in FIG. 1 . Box 
260 reduces fuel data and prepares input for the system simulator. Reduction of 
fuel data involves combining the newly computed fuel with fuels which have 
constant chemistries. Box 270 is system data as indicated; other minor data could 

30 be updated periodically while on-line. Box 280 is the system simulator which, 
given specification of a composite fuel, produces fiiel flow and system efficiency. 
Box 285 is the computation of fuel diemistry and fuel heating value for both the 
composite fuel, and given the fixed compositions of the secondaiy fuels, the 
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composMon of the primary fuel. Box 287 tests for convetgeace of the process 
based on the ccmiposite fuel moles, heating value, and minor t^ms. In general, 
convergences are within 0.5x10"* percent for the computed As-Fired fuel moles. 
Box 288 represents results from the error analysis indicating which instruments 
S require servicing, and numerical minimization of instrummtation otots associated 
with the input measurements of Box 270 versus those computed from Box 280. 
291 repres^ts the transfer of correction factors produced from error analysis for 
effluent instrumentation and/or data assumptions. Box 290 represents a decision 
as to whether the correction factors warrant update and subsequent use. Box 292 

10 determines effluent flow and emission rates, refer to Eqs.(69) thru (71) and 
associated discussion. Box 294 provides the results of the Second Law analysis 
of the thermal ^stem, other ou^ut and reports to system operators such that 
corrective actions rmy take place; also, reports are provided to regulatory 
authorities. Box 296 is a decision to turn the process off (quit) or not. Box 298 

IS is to quit 

FIG. S presents a sample of calcuiational results associated vafh a power 
plant burning coal with wildly varying chemistiy, as discussed earlier in relation 
to FIG. 3. This plant routinely uses coal wiiich varies by ±12 percent in heating 
value. As may be seen in FIG. S» the analysis produced average errors less than 

20 ±0.5 percent in As-Fired heating value, a rraiarkably good result 

FIG. 6 is a schematic representation of a conventional or fluidized bed 
power plant illustrating use of stoichiometric relationships important in applying 
this invention to practical power plants. It should be studied in conjimction with 
Eq.(29) and its associated discussion. Specific discussion of FIG. 6 is presented 

25 earlier. Limestone injection is shown in FIG. 6 which is commonly used in 
fluidized bed power plants. 

In summary, by mass and energy balances based on unity fuel flow rate, 
by usmg accurate molecular weights, by using accurate thermodynamic properties 
of water, steam and combustion gases, by knowing the gross electrical and/or 

30 mechanical power produced from the combustion gases directly or indirectly, by 
knowing the net energy flow supplied to the working fluid, and by recognizing the 
integral relationship of effluent CO2, H2O, SO^, O2, Air/Fuel ratio, air pre-heater 
leakage and O2. concentration in the local combustion air - to their fuel 
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constituents, revised fuel dimustiy masy be computed, ttie heating value based on 
the revised fuel chemistry may be computed, and fuel flow to the system may be 
computed Knowing fuel flow and heating value allows determinadon of effluent 
flows and emission rates. Knowing fuel flow, heating value and fuel chemistry 

S allows routine tracking of a fossil-fired plants* overall thermal efficiency, thus 
continuous correction of problems impacting thermal efficiency is possible. 

Although a particular embodiment of the invention has been described in 
detail for purposes of illustration, various modifications and enhancements may 
be made without departing from the spirit and scope of the invention. 

10 Accordingly, the invention is not to be limited except as by the q)pended claims. 
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CLAIMS 

What is claimed is: 

1 . A method for quantifying the operation of a fossil-fiieled thermial 
system having a heat exchangers/combustion region producing combustion 

5 products, the method comprising the steps of: 
before on-line operation, the steps of 

obtmning reference fuel characteristics, and 
developing explicit mathematical models of the combustion process 
involving at least stoichiometric balances; and thereafter 
10 operating on-line, the step of operating on-line including the steps of 

measuring a set of measurable operating parameters, including at 
least effluent concentrations of 02and CO2, these measurements being made at 
a location downstream of the heat exchangers/combustion region of the thermal 
system, 

IS obtaining an effluent concentration of H2O, if reference fuel 

characteristics indicate fuel water is not predictable, as an obtained effluent H2O, 
obtaining an indicated Air/Fuel ratio, if reference fuel 
characteristics indicate fuel ash is not predictable, as an obtained Air/Fuel ratio, 
obtaining an ambient concentration of O2, 
20 obtaining an air pre-heater leakage factor, and 

calculating a complete As-Fired fuel chemistry, including fiiel 
water and fiiel ash, as a function of the reference fuel characteristics, explicit 
mathematical models of the combustion process, the set of measurable operating 
parameters, the obtained effluent H2O, the obtained Air/Fuel ratio, the ambient 
25 concentration of O2, and the air pre-heater leakage factor. 

2. The method of claim 1, wherein the step of operating on-line 
includes the additional step after calculating the complete As-Fired fuel chemistry, 
of 

calculating consistent moisture-ash-free, dry and As-Fired fuel heating 
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values as a function of the complete As-Fired liiel chemistry and the reference 
fuel characteristics. 

3. The method of claim 2, including, after the step of calculating 
consistent moisture-ash-free, dry and As-Fired heating values, the additional steps 

5 of 

obtaining fuel factors associated with the As-Fired fuel chemistry, 
^ updating the reference fuel characteristics based on the calculated As-Fired 
fuel chemistry and its sensitivities to fuel characteristics, 

completing a multidimensional minimization analysis employing fiiel 
1 0 factors and heating values to minimize the collective error associated with at least 
one of the measured effluent CO2, the obtained effluent H2O, the indicated 
Air/Fuel ratio, the ambient concentration of O2 and the air pre-heater leakage 
factor, 

calculating and applying for subsequent on-line analysis correction factors 
15 to the measured effluent CO2, obtained efiQuent H2O, the indicated Air/Fuel ratio, 
the ambient concentration of O2, and the air pre-heater leakage factor. 

4. The method of claim 2, vrtierein the set of measurable operating 
parameters includes effluent temperature, and wherein the method includes an 
additional step, after the step of calculating consistent As-Fired fiiel heating value, 

20 of 

calculating a consistent boiler efficiency as a function of the complete As- 
Fired fuel chemistiy, efflumt temperature, effluent concentrations, and heating 
value. 

5. The method of claim 4, wherein the step of operating on-line 
25 includes an additional step of 

measuring the useful output derived directly from the combustion products, 
and the step of measuring a set of measurable operating parameters includes the 
step of 

measuring a useful power and energy flow produced from a working fluid 



CA. 02325929 2003-11-20 



-65- 

heated by combustion products, 

determining an efficiency associated with the working fluid's process, and 
further including an additional step, aft^ the step of calculating a consistent boiler 
efficiency, of 

5 calculating a system thermal efficiency as a function of the efficiency 

associated with the working fluid's process and the consistent boiler efficiency. 

6. The method of claim 4, wherein the step of operating on-line 
includes an additional step of 

measuring a useflil output derived directly from the combustion products, 
10 and including an additional step, after the step of calculating the complete As- 
Fired fuel chemistry, of 

calculating a consistent heating value as a function of the complete As-Fired 
fuel chemistry, and 

including an additional step, after the step of calculating a consistent heating value, 
15 of 

calculating a consistent As-Fired fiiel flow rate as a function of the 
consistent boiler efficiency, the useful ou^ut, and the As*Fired heating value. 

7. The method of claim 6, including an additional step, after the step 
of calculating a consistent As-Fired fuel flow rate and boiler efficiency, of 

20 calculating a consistent total effluent flow rate as a fimction of the 

consistent As-Fired fiiel flow rate. 

8. The method of claim 7, including an additional step, after the step 
of calculating a consistent total effluent flow rate, of 

obtaining a set of effluent concentrations including at least SO2 and NOx» 
25 calculating individual effluent flow rates as fiinctions of the total effluent 

flow rate and effluent concentrations. 

9. The method of claim 8, including an additional step, after the step 
of calculating individual effluent flow rates, of 
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calculating consistent individual emission rates as a function of the 
individual efiSuent flow rates, consistent As-Fired fuel flow rate, and consistent As- 
Fired heating value. 

10. The method of claim 1, wherein the thermal system comprises a 
5 thermal system selected from the group consisting of a coal-buming power plants, 
oil-burning power plants, gas-fired power plants, biomass combustors, fluidized 
bed combustors, a conventional electric power plant, a steam generator, a package 
boiler, a combustion turbine, and a combustion turbine with a heat recovery boiler. 

1 L The method of claim 8, including an additional step, after the step 
10 of calculating consistent individual effluent flow rates and those quantities on 
which they are consistency based, of 

calculating a Second Law analysis of the thermal system and Fuel 
Consumption Indices. 

12. A method for quantifying the operation of a fossil-fiieled thermal 
1 S system having a heat exchangers/combustion region wherein a fuel flow produces 

combustion products, the method comprising the steps of: 

operating on-line, the step of operating on-line including the steps of 

obtaining an As-Fired fuel chemistry including at least one fuel 
constituent selected from the group consisting of weight percent carbon, weight 
20 percent hydrogen, and weight percent oxygen of the fiiel; and 

calculating a heating value of the fiiel in Btu per pound of fiiel, the 
step of calculating including a step of forming products of numerical coefficients 
times the weight percent of the fiiel constituent, wherein at least one of the 
numerical coefficients is selected from the group consisting of from 145 to 1 54 for 
25 weight percent carbon, from 385 to 570 for weight percent hydrogen, and from -5 1 
to -57 for weight percent oxygen, and combinations thereof. 

13. A method for quantifying the operation of a fossil-fiieled fluidized 
bed thermal system producing combustion products, the method comprising the 
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steps of: 

before on-line operation, the steps of 

obtaining reference fuel characteristics, and 
developing explicit mathematical models of the combustion process 
S involving at least stoichiometric balances; and thereafter 

operating on-line, the step of operating on-line including the steps of 

measuring the measurable operating parameters, including at least 
an effluent temperature, and 

the concentration of effluent O2, CO2, and SO2, at a location 
10 downstream of a heat exchangers/combustion region of the fluidized bed thermal 
system; 

measuring an indicated limestone flow rate; 

measuring a useful ou^ut derived directly from the combustion 

products, 

15 calculating a complete As-Fired fuel chemistry, wherein the 

complete As-Fired fuel chemistry includes the elementary constituents, fuel water, 
and fiiel ash of the ftiel, and 

determining the Eneigy Released during Combustion based on 
explicit mathematical models of the combustion process and the complete As-Fired 
20 fuel chemistry, 

determining a consistent fuel flow, a consistent effluent flow, and 
a consistently calculated system efficiency, as a function of the measurable 
operating parameters, the indicated limestone flow rate, the useful output, and the 
complete As-Fbed fuel chemistry. 

25 14. A method for quantifying the operation of a fossil- fueled thermal 

system in which a fossil fiiel is supplied at a flow rate to a heat 
exchangers/combustion region and combusted to produce hot combustion gases, 
which heats a working fluid and/or produces shaft power then exits through an 
exhaust stack, the method comprising the following steps: 

30 performing an off-line operation comprising the steps of 

obtaining reference fuel characteristics. 
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obtaining reference fuel stoichiometric data, 

obtaining current measurements of the system*s operating 
parameters, and 

performing an on-line operation comprising the steps of 

measuring the useful output of the system, 

obtaining fuel data and characteristics, the step of obtaining fuel data 
including the step of obtaining composite fuel concentrations and composite 
heating value, if multiple fuels are used, 

introducing fuel concentrations and heating values to a mathematical 
model of the thermal system, 

obtaining routine systems operational parameters, 

obtaining values of the effluents O2, CO2, CO, H2O, SO2 and NOx, 

obtaining the ambient concentration of O2, 

obtaining air pre-heater leakage and dilution factors, 

obtaining an indicated Air/Fuel ratio, 

computing molar moisture-ash-firee fiactions of fuel carbon, fiiel 
water and fuel ash as explicit stoichiometric solutions, dependent at least in part on 
the effluents O2, CO2, SO2, H2O, Air/Fuel ratio, air pre-heater factors and ambient 
concentration of 62, 

finding the molar moisture-ash-fi*ee fractions of fuel nitrogen, 
oxygen, hydrogen, and sulfur, 

converting the molar moisture-ash-fi:ee fiiel concentrations to a 
molar dry base, then to a molar As-Fired wet base, and finally to As-Fired wet 
weight fractions, to obtain a complete and consistent computed As-Fired fuel 
chemistry, 

computing a heating value based on a moisture-ash-fiee weight base, 
then converted to a diy base, and th^ to a weight-based As-Fired heating value, 
and 

executing the mathematical model of the thermal system using the 
fiiel information and the concentration of effluent O2 to produce consistent 
stoichiometric values of effluent COj, SO2 and H2O values, the Air/Fuel ratio, the 
moles of fiiel per basis moles of dry gaseous effluent, and at least the following 
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self-consistent thermal performance parameters: As-Fired fuel flow, effluent flow, 
emission rates, boiler efficiency, and over-all system thermal efficiency. 

15. The method of claim 14, including an additional step, after the step 
of executing, of 

performing analysis of instrumentation errors to obtain correction factors, 
and, if excessive, applying the correction factors to instrumentation signals such 
that subsequent on-line operation produces minimum errors in fuel chemistry and 
heating value determinations. 

16. The method of cla'mi 13, including an additional step, after the step 
of determimng a consistent fiiel flow, of 

determining a consistent effluent flow, and 
reporting the consistent effluent flow. 

17. The method of claim 13, including an additional step, after die step 
of determining a consistent fiiel flow, of 

calculating an effluent, ash flow, and 
reporting the effluent ash flow. 

1 8. The method of claim 13, including an additional step, after the step 
of determining a consistent fuel flow, of 

calculating an set of emission rates, and 
reportmg the set of emission rates. 

19. The method of claim 13, including an additional step, after the step 
of determining an Energy Released during Combustion, of 

reporting to the operator of the fossil-fueled fluidized bed thermal system 
the Energy Released during Combustion. 

20. A method for quantifying the operation of a fossil-fueled thermal 
system having a heat exchangers/combustion region producing combustion 
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products, the method comprising the steps of: 
before on-line operation, the steps of 

obtaining a set of reference fuel characteristics, and 

developing explicit mathematical models of the combustion process 

involving at least stoichiometric balances; and thereafter 

operating on-line, the step of operating on-line including the steps of 

measuring a set of measurable operating parameters, including at 

least effluent concentrations of O2 and CO2, these measurements being made at a 

location downstream of the heat exchangers/combustion region of the thermal 

system, 

obtaining an effluent concentration of H2O if the set of reference fuel 
characteristics indicates that fuel water is not predictable, as an obtained effluent 
H2O, 

obtaining a concentration of O2 in the ambient air entering the 
thermal system, 

obtaining an air pre-heater leakage factor, 

calculating a set of fiiel chemistry concentrations including 
elementary fiiel constituents and fuel water, but excluding fuel ash, as a function 
of the set of reference fuel characteristics, explicit mathematical models of the 
combustion process, the set of measurable operating parameters, the obtained 
efnuent HjO, the concentration of O2 in the ambient air entering the thermal 
system, and the air pre-heater leakage factor. 

21. The method of claim 20, wherein the step of operating on-line 
includes the additional steps after the step of calculating the set of fiiel chemistry 
concentrations including elementary fiiel constituents and fuel water, of 

obtaining an indicated fiiel flow if the set of reference fiiel 
characteristics indicates that fuel ash is not predictable, as an obtmned fuel flow 
used for calculating fuel ash concentration, 

calculating a fiiel ash concentration if the set of reference fuel 
characteristics indicates fuel ash is not predictable, as a function of the set of fuel 
chemistry concentrations including elementary fuel constituents and fuel water, 
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. explicit mathematical models of die combustion process, the obtained fuel flow, 
the concentration of O2 in the ambient air entering the thermal system, and the air 
pre*heater leakage factor, and 

calculating a complete As-Fired fiiel chemistty as a function of the 
set of fuel chemistiy concentrations which includes elementaiy fuel constituents, 
fuel water and fuel ash concentrations. 

22. The method of claim 20, wherein the step of operating on-line 
includes the additional steps after the step of calculating the set of fuel chemistry 
concentrations including elementary fiiel constituents and fuel water, of 

obtaining a normalized Air/Fuel ratio if the set of reference fuel 
characteristics indicates that fiiel ash is not predictable, as an obtained Air/Fuel 
ratio used for calculating fuel ash concentration, 

calculating a fuel ash concentration if the set of reference fiiel 
characteristics indicates fiiel ash is not predictable, as a function of the set of fiiel 
chemistry concentrations including elementary fuel constituents and fiiel water, 
explicit mathematical models of the combustion process, the obtained Air/Fuel 
ratio, the concentration of O2 in the ambient air entering the thermal system, and 
the air pre-heater leakage factor, and 

calculating a complete As-Fired fiiel chemistiy as a function of the 
set of fiiel chemistiy concentrations which includes elementary fuel constituents, 
fiiel water and fiiel ash concentrations. 

23. The method of claim 20, wherein the step of operating on-line 
includes the additional steps after the step of calculating the set of fiiel chemistry 
concentrations including elementary fuel constituents and fiiel water, of 

obtaiiiing a predictable fiiel ash concentration based on the set of 
reference fuel characteristics, and 

calculating a complete As-Fired fiiel chemistry as a fimction of the 
set of fiiel chemistry concentrations which includes elementary fiiel constituents, 
fiiel water and fiiel ash concentrations. 
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24. The method of claim 20, wherein the step of calculating the set of 
fuel diemistry concentrations including elementary fuel constituents and fuel 
water, includes the step of 

calculating a set of moisture-ash-free fuel chemistry concentrations 
5 including elementary fuel constituents and fiiel water, but excluding fuel ash, as a 
function of the set of reference fuel characteristics, explicit mathematical models 
of the combustion process, the set of measurable operating parameters, the 
obtained effluent HjO, the concentration of O2 in the ambient air entering the 
thermal system, and the air pre*heater leakage factor. 

10 25. The method of claim 20, >\1ierein the step of calculating the set of 

fuel chemistry concentrations including elementary fuel constituents and fuel 
water, mcludes the step of 

calculating a set of dry-based fiiel chemistry concentrations 
including elementary fuel constituents and fiiel water, but excluding fiiel ash, as a 

1 S function of the set of reference fuel characteristics, explicit mathematical models 
of the combustion process, the set of measurable operating parameters, the 
obtained effluent HjO, the concentration of O2 in the ambient air entering the 
thermal system, and the air pre-heater leakage factor. 

26. The method of claim 2 1 , where'm the step of calculating the fiiel ash 
20 concentration, includes the step of 

calculating explicitly a moisture-ash-£ree fiiel ash concentration if 
the set of reference fuel characteristics indicates fiiel ash is not predictable, as a 
fiinction of the set of fiiel chemistiy concentrations including elementary fiiel 
constituents and fiiei water, explicit mathematical models of the combustion 
25 process, the obtained fiiel flow, the concentration of O2 in the ambient air entering 
the thermal system, and the air pre-heater leakage factor. 



27. The method of claim 2 1 , wherein the step of calculating the fuel ash 
concentration, includes the step of 

calculating explicitly a dry-based fuel ash concentration if the set of 



\ 
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reference fuel characteristics indicates fuel ash is not predictable, as a function of 
the set of fiiel chemistty concentrations including elementary fiiel constituents and 
fuel water, explicit mathematical models of the combustion process, the obtained 
fiiel flow, the concentration of O2 in the ambient air entering the thermal system, 
5 and the air pre-heater leakage factor. 

28. The method of claim 22, wherein the step of calculating the fuel ash 
concentration, includes the step of 

calculating explicitly a moisture-ash-free fuel ash concentration if 
the set of reference fuel characteristics indicates fiiel ash is not predictable, as a 
10 fimction of the set of fiiel chemistry concentrations including elementary fuel 
constituents and fiiel water, explicit mathematical models of the combustion 
process, the obtained Air/Flow ratio, the concentration of O2 in the ambient air 
entering the thermal system, and the mr pre-heater leakage factor. 

29. The method of claim 22, wh^in the step of calculating the fuel ash 
1 S concmtration, includes the step of 

calculating explicitly a diy-based fiiel ash concentration if the set of 
reference fuel characteristics indicates fiiel ash is not predictable, as a fimction of 
the set of fiiel chemistry concentrations including elementary fiiel constituents and 
fiiel water, explicit mathematical models of the combustion process, the obtained 
20 Air/Fuel ratio, the concentration of O2 in the ambient sdr entering the thermal 
system, and the air pre-heater leakage factor. 

30. The method of claim 23, wherein the step of obtmning the 
predictable fiiel ash concentration based on the set of reference fiiel characteristics, 
includes the step of 

25 obtaining a predictable moisture-ash-free fiiel ash concentration 

based on the set of reference fiiel characteristics. 

31. The method of claim 23, wherein the step of obtaining the 
predictable fiiel ash concentration based on the set of reference fiiel characteristics. 
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includes the step of 

obtaining a predictable diy-based fuel ash concentration based on 
the set of reference fiiel characteristics. 

32. The method of claim 21, wherein the step of operating on-line 
5 includes the additional step after calculating the complete As-Fired fiiel chemistiy, 

of 

calculating an As-Fired fiiel heating value as a fimction of the complete As- 
Fired fiiel chemistry and the set of reference fiiel characteristics. 

33. The method of claim 32, including, after the step of calculating the 
1 0 As-Fired ftiel heating value, the additional steps of 

obtaining fiiel factors associated with the complete As-Fired fiiel chemistry, 
completing a multidimensional minimization analysis employing fiiel 
factors and heating values to minimize die collective error associated with at least 
one of the measured effluent CO2, the obtained effluent H2O, the obtained fiiel 
1 5 flow, the concentration of O2 in the ambient air entering the thermal system, and 
the air pre-heater leakage factor, 

obtaining and applying for subsequent on-line analysis correction factors 
to the measured effluent CO2, the obtained effluent H2O, the obtained fiiel flow, 
the concentration of O2 in the ambient air entering the thermal system, and the air 
20 pre-heater leakage factor. 

34. The method of claim 32, wherein the set of measurable operating 
parameters includes effluent temperature, and wherein the method includes an 
additional step, after the step of calculating the As-Fired fiiel heating value, of 

obtaining a system energy credit term, 
25 calculating a consistent boiler efficiency as a fimction of the complete As- 

Fired fiiel chemistry, effluent temperature, the effluent concentrations, the As-Fired 
fiiel heating value and the system energy credit term. 
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35. The method of claim 34, wherein the step of operating on-line 
includes an additional step of 

measuring the useful ou^ut derived directly from the combustion products, 
and the step of measuring a set of measurable operating parameters includes the 
5 step of 

measuring a usefiil electrical and thermal power produced from a woricing 
fluid heated by combustion products, 

determining an efiFiciency associated with the working fluid's process, and 
further including an additional step, after the step of calculating a consistent boiler 
10 efficiency, of 

calculating a system thermal efficiency as a function of the efficiency 
associated with the working fluid's process and the consistent boiler efficiency. 

36. The method of claim 34, wherein the set of measurable operating 
parameters includes measuring a useful output derived directly from the 

1 5 combustion products, and including an additional step, after the step of calculating 
the consistent boiler efficiency, of 

calculating a consistent As-Fired fuel flow rate as a function of the 
consistent boiler efficiency, the useful output, the As-Fired fuel heating value and 
the system energy credit term. 

20 37. The method of claim 36, including an additional step, after the steps 

of calculating the consistent As-Fired fuel flow rate and the consistent boiler 
efficiency, of 

calculating a consistent effluent flow rate as a function of the consistent As- 
Fired fiiel flow rate. 

25 38. The method of claim 37, including an additional step, after the step 

of calculating the consistent eSluent flow rate, of 

obtaining a set of effiuent concentrations including at least SO2 and NOx, 
calculating a set of individual effluent flow rates including those of at least 

SO2 and NOx as functions of the total effiuent flow rate and the set of effiuent 
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concentrations. 

39. The method of claim 38, including an additional step, after the step 
of calculating the set of individual effluent flow rates, of 

calculating a set of consistent emission rates as a function of the individual 
S effluent flow rates, consistent As-Fired fiiel flow rate and consistent As-Fired 
heating value. 

40. TTie method of claim 21, wherein the thermal system comprises a 
thermal system selected from the group consisting of a coal-burning power plants, 
oil-burning power plants, gas-fired power plants, biomass combustors, fluidized 

1 0 bed combustors, a conventional electric power plant, a steam generator, a package 
boiler, a combustion turbine, and a combustion turbine with a heat recovery boiler. 

4 1 . The method of claim 37, including an additional step, after the step 
of calculating the consistent effluent flow rate , of 

calculating a Second Law analysis of the thermal system and Fuel 
IS Consumption Indices. 

42. The method of claim 22, wherein the step of operating on-line 
includes the additional step after calculating the complete As-Fired fiiel chemistry, 
of 

calculating an As-Fired fuel heating value as a function of the complete As- 
20 Fired fuel chemistry and the set of reference fiiel characteristics. 

43. The method of claim 42, including, after the step of calculating the 
As-Fired fiiel heating value, the additional steps of 

obtaining fuel factors associated with the complete As-Fired fuel chemistry, 
completing a multidimensional minimization analysis employing fuel 
25 factors and heating values to minimize the collective error associated with at least 
one of the measured effluent CO2, the obtained eflnuent H2O, the obtained fuel 
flow, the concentration of O2 in the ambient air entering the thermal system, and 
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the air pre-heater leakage factor, 

obtaining and applying for subsequent on-line analysis correction factors 
to the measured efiQuent CX)2, the obtained effluent H2O, the obtained fuel flow, 
the concentration of O2 in the ambient air entering the thermal system, and the air 
5 pre-heater leakage factor 

44. The method of claim 42, wherein the set of measurable operating 
parameters includes effluent temperature, and wherein the method includes an 
additional step, after the stq) of calculating the As-Fired fiiel heating value, of 

obtaining a system energy credit term, 
1 0 calculating a consistent boiler efGciency as a function of the complete As- 

Fired fuel chemistry, effluent temperature, the effluent concentrations, the As-Fired 
fuel heating value and the system energy credit term. 

45. The method of claim 44, wherein the step of operating on-line 
includes an additional step of 

1 S measuring tiie useful ou^ut derived (tirectly firom the combustion products, 

and the step of measuring a set of measurable operating parameters includes the 
step of 

measuring a useful electrical and thermal power produced from a working 
fluid heated by combustion products, 
20 determining an efficiency associated with the woridng fluid's process, and 

further including an additional step, after the step of calculating a consistent boiler 
efflciency, of 

calculatii^ a system thermal efficiency as a function of the efficiency 
associated with the working fluid's process and the consistent boiler efficiency. 

25 46. The method of claim 44, where'm the set of measurable operating 

parameters includes measuring a useful output derived directly from the 
combustion products, and including an additional step, after the step of calculating 
the consistent boiler efficiency, of 

calculating a consistent As-Fired fuel flow rate as a function of the 
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consistent boiler efficiency, the useful output, the As-Fired fuel heating value and 
the system energy credit term. 

47. The method of claim 46, including an additional step, after the steps 
of calculating the consistent As-Fired fuel flow rate and the consistent boiler 

5 efficiency, of 

calculating a consistent effluent flow rate as a function of the consistent As- 
Fired fiiel flow rate. 

48. TTie method of claim 47, including an additional step, after the step 
of calculating the consistent effluent flow rate, of 

1 0 obtaining a set of effluent concentrations including at least SO2 and NOx, 

calculating a set of individual effluent flow rates including those of at least 
SO2 and NOx ^ functions of the total effluent flow rate and the set of effluent 
concentrations. 

49. The method of claim 48, including an additional step, after the step 
1 S of calculating the set of individual effluent flow rates, of 

calculating a set of consistent emission rates as a function of the individual 
effluent flow rates, consistent As-Fired fiiel flow rate and consistent As-Fired 
heating value. 

50. The method of claim 22, wherein the thermal system comprises a 
20 thermal system selected fiom the group consisting of a coal-buming power plants, 

oil-burning power plants, gas-fired power plants, biomass combustors, fluidized 
bed combustors, a conventional electric power plant, a steam generator, a package 
boiler, a combustion turbine, and a combustion turbine with a heat recovery boiler. 

5 1 . The method of claim 47, including an additional step, after the step 
25 of calculating the consistent effluent flow rate , of 

calculating a Second Law analysis of the thermal system and Fuel 
Consumption Indices. 
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52. The method of claim 23, wherein the step of operating on-line 
includes the additional step after calculating the complete As-Fired fiiel chemistry^ 
of 

calculating an As-Fired fuel heating value as a function of the complete As- 
S Fired fuel chemistry and the set of reference fuel characteristics. 

53. The method of claim 52, including, after the step of calculating the 
As-Fired fiiel heating value, the additional steps of 

obtaining fuel factors associated with the complete As-Fired fiiel chemistry, 
completing a multidimensional minimization analysis employing fuel 
10 factors and heating values to minimize the collective error associated with at least 

one of the measured effluent CO2, the obtained effluent H2O, the obtained fuel 

flow, the concentration of O2 in the ambient air entering the thermal system, and 

the air pre-heater leakage factor, 

obtaining and applying for subsequent on-line analysis correction factors 
15 to the measured effluent CO2, the obtained effluent H2O, the obtained fuel flow, 

the concentration of O2 in the ambient m entering the thermal system, and the mr 

pre-heater leakage factor. 

54. The method of claim 52, wherein the set of measurable operating 
parameters includes effluent temperature, and wherein the method includes an 

20 additional step, after the step of calculating the As-Fired fiiel heating value, of 
obtaining a system energy credit term, 

calculating a consistent boiler efficiency as a function of the complete As- 
Fired fiiel chemistty , effluent temperature, the effluent concentrations, the As-Fired 
fiiel heating value and the system energy credit term. 

25 55. The method of claim 54, wherein the step of operating on-line 

includes an additional step of 

measuring the useful output derived directly from the combustion products, 
and the step of measuring a set of measurable operating parameters includes the 
step of 
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measuring a useful electrical and thermal power produced from a working 
fluid heated by combustion products, 

determining an efficiency associated with the working fluid's process, and 
further including an additional step, after the step of calculating a consistent boiler 
5 efficiency, of 

calculating a system thermal efficiency as a function of the efficiency 
associated with the working fluid's process and the consistent boiler efficiency. 

56. The method of claim 54, wherein the set of measurable operating 
parameters includes measuring a useful output derived directly from the 
1 0 combustion products, and including an additional step, after the step of calculating 
. the consistent boiler efficiency, of 

calculating a consistent As-Fired fuel flow rate as a function of the 
consistent boiler efficiency, the useful output, the As-Fired fuel heating value and 
the system energy credit term. 

1 5 57. The method of claim 56, including an additional step, after the steps 

of calculating the consistent As-Fired fuel flow rate and the consistent boiler 
efficiency, of 

calculating a consistent effluent flow rate as a fimction of the consistent As- 
Fired fuel flow rate. 

20 58. The method of claim 57, including an additional step, after the step 

of calculating the consistent effluent flow rate, of 

obtaining a set of effluent concentrations including at least SO2 and NOx, 
calculating a set of individual effluent flow rates including those of at least 
SO2 and NOx ^ functions of the total effluent flow rate and the set of effluent 
25 concentrations. 



59. The method of claim 58, including an additional step, after the step 
of calculating the set of individual effluent flow rates, of 

calculating a set of consistent emission rates as a function of the individual 
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effluent flow rates, consistent As-Fired fuel flow rate and consistent As-Fired 
heating value. 

60. The method of claim 23, wherein the thermal system comprises a 
thermal system selected from the group consisting of a coal-burning power plants, 

5 oil-burning power plants, gas-fired power plants, biomass combustors, fluidized 
bed combustors, a conventional electric power plant, a steam generator, a package 
boiler, a combustion turbine, and a combustion turbine with a heat recovery boiler 

6 1 . The method of claim 5 7, including an additional step, after the step 
of calculating the consistent effluent flow rate , of 

10 calculating a Second Law analysis of the thermal system and Fuel 

Consumption Indices. 

62. The method of claim 1 , wherein the step of calculatmg the complete 
As-Fired fuel chemistry includes the step of 

calculating explicitly a complete As-Fired fiiel chemistry, including fuel 
IS water and fuel ash, as a function of the reference fuel characteristics, explicit 
mathematical models of the combustion process, the set of measurable operating 
parameters, the obtained effluent H2O, the obtained Air/Fuel ratio, the 
concentration of O2 in the ambient air entering the thermal system, and the air pre- 
heater leakage factor. 

20 63 . The method of claim 1 , wherein the step of obtiuning the indicated 

Air/Fuel ratio, if reference fiiel characteristics indicate fuel ash is not predictable, 
as the obtained Air/Fuel ratio, includes the step of 

obtaining a normalized Air/Fuel ratio, if reference fuel characteristics 
indicate fuel ash is not predictable, as an obta'med Air/Fuel ratio. 

25 64. The method of claim 1 , wherein the step of obtaining the ambient 

concentration of O2 includes the step of 

using a value of 20.948 percent for the ambient concentration of oxygen. 
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65. The method of claim 1 , wherein the step of obtaining the ambient 
concentration of oxygen includes the step of 

using an average value at sea level determined by the Unites States' 
National Aeronautics and Space Administration for the ambient concentration of 
5 oxygen. 

66. Hie method of claim 1, wherein the step of obtaining the air pre- 
heater leakage factor includes the step of 

using a value of unity for the air pre-heater leakage factor. 

67. The method of claim 20, wherein the step of calculating the set of 
10 fuel chemistry concentrations including elementary fuel constituents and fuel 

water, but excluding fiiel ash, includes the step of 

calculating explicitiy aset of fuel chemistry concentrations including 
elementaiy fuel constituents and fuel water, but excluding fuel ash, as a fimction 
of the set of reference fuel characteristics, explicit mathematical models of the 
15 combustion process, the set of measurable operating parameters, the obtained 
effluent H2O, the concentration of O2 in the ambient air entering the thermal 
system, and the air pre-heater leakage factor. 

68. The method of claim 20, wherein the step of obtaining the 
concentration of O2 in the ambient air entering the thermal system includes the step 

20 of 

using a value of 20.948 percent for the concentration of O2 in the ambient 
mr entering the thermal system. 

69. The method of claim 20, wherein the step of obtaining the 
concentration of O2 in the ambient air entering the thermal system includes the step 

25 of 

using an average value at sea level determined by the United States' 
National Aeronautics and Space Administration for the concentration of O2 in the 
ambient air entering the thermal system. 
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70. The method of claim 20, wherein the step of obtaining the air pre- 
heater leakage factor includes the step of 

using a value of unity for the air pre-heater leakage factor. 
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